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(54) Tide: METHOD FOR CONTROLLING CHEMICAL SYNTHESIS PROCESSES 

(54) Utre: METHODE DE REGULATION DE PROCEDES DE SYNTHESE DE PRODUITS CHIMIQUES 
(57) Abstract 

A method for controlling a process for synthesising at least one chemical in 
a plant including at least one reactor (R) that may be considered as a perfectly 
mixed reactor, wherein control variables (GC) enable the process to be controlled 
so that a number of variables known as controlled variables (GR) and related to the 
properties of the product and/or the operation of the process match corresponding 
set values (Cgr). Said method comprises the steps of (a) inputting set values for the 
controlled variables (Cgr); (b) using a prediction unit (OP) to calculate predictions 
of the controlled variables (Pgr) on the basis of measurements of the process control 
vLSKS (O using a contiol unit (OC) to calculate set values for the process control variables (Codon the basis of the set values 
(^)and pSdictions (Pgr) of the controlled variables; and (d) transmitting set values for the process contro variables (CccMo process 
coSrol actuators. The prediction unit (OP) is based on a mathematical process model known as a direct model (M). and designed m such 
a way that the mass Mxr of at least one component (X) in the reactor (R) is predicted by means of the equanon Mxr " ^ffxato . rx. 
rJ) wherein Fxwn is the mass flow rate of component X at the inlet of the reactor R. r x is the dwell time of x m the reartor, and the 
fZiZ7y-L%?£ r) is the solution of the differential equation r . dy/dt + y. calculated with the instantaneous value of u and r as 
well as with die last calculated value of y. 

(57) Abrlgtf 

Une mtfthode de regulation d'un proctfdtf de synthese d'au moins un produit chimique dans un appareillage comprer^t au moins un 
reacteur (R) pouvant etre assimiltf a un reacteur parfaitement melange dans laquelle des grandeurs de comrnande (GC g"^' d >girtur 
kderoulement du precede en vue que des ou plusieurs grandeurs liees aux propritftes du produit et/ou au denotement <hj precede, appelfes 
ZSn^Sta (GR). soient tfgales a des consign* coirespondantes (Cgr). ladite ^^^^rS £Ld e * 
dVconsignes concernant les grandeurs regltfes (Cgr); (b) calcul. au moyen d'un organe de prediction (OP), de predicUonsdes l£nd«» 
rfritfes (Erf sur la base de mesures des grandeurs de comrnande du procedtf (Mac); (c) utilisation d'un organe de contrtle (OQ pour 
calcul^r de?~^gL des grandeurs de comrnande du precede (Cgc). sur base des consignes (Cor) et des predictions (Pgr) des g^dcura 
W uSTssion defconsignes des grandeurs de comrnande du procedtf (Cgc) a des actionneurs afin d'agn ^sur le^roulement du 
prS; dans laquelle 1'organe de prediction (OP) est bastf sur un modtfle mathtfmatique du proctfdtf, appeltf modele direct <M), et est concu 
SmTftreVon prtfdit la masse Mxr d'au moins un constimant (X) dam le reacteur (R) par 1 'equanon: Mxr |J^M^r-^ rx. 
^ dans laouelle Fxrun est le dtfbit massique du constituant X entrant dans le reacteur R; rx est le temps de stfjour de x dans le reacteur 
[alonc^nT^llGfe r) est la Solution de 1'tfquation difftfrentielle r . dy/dt + y. calculee avec la valeur instantantfe de u et de r. ains. 
qu'avec la demiere valeur de y calculeee. 
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Methode de regulation de precedes de synthese de produits chimicmes 



La presente invention concerne une methode de regulation de 
precedes de synthese de produits chimiques. EUe concerne egalement un 
dispositif de regulation pour la mise en oeuvre de cette methode, ainsi qu'un 
precede de synthese, en particulier de polymere, regule par cette methode, 
5 Dans un procede de synthese de produits chimiques conduit 

conventionnellement, des regulateurs de type PID (proportionnel - integral - 
differentiel) sont utilises pour reguler individuellement un nombre plus ou 
moins important de grandeurs (temperatures, debits, pressions, ...) ayant 
une influence sur le derouiement de la synthese. En d'autres termes, pour 

10 chaque temperature, debit ou pression a reguler, on mesure en continu (ou 
par intermittence) sa valeur effective, et un regulateur PID compare cette 
valeur effective a une valeur de consigne et agit sur la grandeur a 
commander de faqon a reduire, le cas echeant, la difference entne la valeur 
de consigne et la valeur mesuree. 

15 Vu la complexity de la plupart des precedes industriels de synthese 

chimique, les valeurs de consignes des differents regulateurs doivent 
aujourd'hui encore etre ajustees empiriquement pour obtenir finalement les 
proprietes souhaitees du produit synthetise. On utilise a cette fin des 
recettes, qui fournissent des combinaisons de parametres determiners 

20 empiriquement pour obtenir, en regime stable, les proprietes souhaitees du 
produit synthdtise. 

De ces recettes pourraient etre deduites, k Taide d'outils statistiques 
plus ou moins sophistiques, des relations empiriques entre les grandeurs 
reglees et les proprietes du produit synthetise. II est cependant evident que 

25 ces relations empiriques ne peuvent guere tenir compte des multiples 

interdependances existant entre les differentes grandeurs reglees separdment, 
pas plus que des perturbations inconnues telles que les teneurs en impuretes 
des matieres premieres. 

II est aussi evident qu'une regulation classique a boucle fermee, uti- 

30 lisant comme corrections de retour des mesures de proprietes essentielles du 
produit synthetise, est difficilement applicable a la plupart des procedes de 
synthese. En effet, les temps morts intervenant soit dans le procede, soit 
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dans les mesures ou analyses utilisees comme corrections de retour, sont 
trop eleves, et les interdependances entre les differentes grandeurs regissant 
le procede sont trop complexes. 

II est connu depuis longtemps que les procedes de synthese, 
5 notamment les procedes de synthase en continu de polymeres (procedes de 
polymerisation), equipes de regulations avec ajustage empirique des valeurs 
de consignes, ont des inconvenients majeurs, qui peuvent etre resumes 
comme suit : 

- les demarrages du procedd de synthese prennent beaucoup de temps et 
10 genferent des quantites importantes de produit "hors normes" ; 

- les transitions de qualite sont lentes, ce qui entraine egalement la produc- 
tion de quantites importantes de produits de transition "hors normes" ; 

- failure de marche du procede, c'est-a-dire la masse de produit(s) 
synth&isee par unite de temps, peut difficilement etre modifiee sans 

15 alterer les proprietes de ce ou ces produits ; 

- la Constance des proprietes essentielles du ou des produits synthetises 
laisse souvent & desirer, meme en regime stable. 

Afin d'eviter un ajustage empirique des valeurs de consigne, il a ete 
propose dans la litterature spec iali see d'utiliser des methodes de regulation 

20 de procedes de synthese qui ont recours a des equations caracteristiques 

modelisant le procede de synthese pour relier certaines proprietes du ou des 
produits synth6tises aux conditions de marche du ou des reacteurs lors de la 
synthese. Cependant, afin de limiter la complexity de ces equations caracte- 
ristiques, on considerait jusqu'a present qu'on etait oblige en pratique soit 

25 de considerer exclusivement le cas statique (regime stable), soit de se 

limiter a une modelisation empirique tres simplifiee de la dynamique du 
procede. L'utilisation d'un modele statique est limitee au controle d'un 
regime de production assez stable. 

Dans le cas d'une modelisation empirique, les equations 

30 caracteristiques ne sont valables que pour une plage de validite etroite (a 
proximite du point oil la modelisation a ete effectuee). Dans les deux cas, 
les phases de demarrage et les phases de transition sont mal maitrisees. II 
est certes envisageable de "couvrir" une plage plus large de conditions de 
fonctionnement en procedant a une telle modelisation locale en plusieurs 

35 points distincts de l'espace des parametres de fonctionnement, mais une telle 
approche devient prohibitive des que Ton tente de reguler plusieurs gran- 
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deurs en agissant sur de multiples parametres. 

II serait par consequent souhaitable de disposer d'une methode et 
d'un dispositif de regulation simples, mieux adaptes aux specificites de la 
dynamique de precedes de synthese de produits chimiques. 
5 A cette fin, la pr^sente invention conceme une methode de regulation 

d'un precede de synthase d'au moins un produit chimique dans un 
appareiilage comprenant au moins un reacteur (R) pouvant etre assimil£ & 
un reacteur parfaitement melange, dans laquelle une ou plusieurs grandeurs 
de commande (GC) permettent d'agir sur le deroulement du precede en vue 
10 qu'une ou plusieurs grandeurs \\6es aux proprietes du produit et/ou au 

deroulement du precede, appelees grandeurs reglees (GR), soient 6gales k 
des consignes conrespondantes (Cq R ) (ou du moins aussi proches que 
possible de celles-ci), ladite methode comprenant les etapes suivantes : 
(a) entree de consignes concernant les grandeurs reglees (C GR ); 
15 (b) calcul, au moyen d'un organe de prediction (OP), de predictions des 
grandeurs reglees (Pgr)* sur ,a base de mesures des grandeurs de 
commande du precede (Mq^); 

(c) utilisation d'un organe de controle (OC) pour calculer des consignes 
des grandeurs de commande du precede (Cqq), sur base des 

20 consignes (Cq R ) et des predictions (Pgr) des grandeurs regimes; 

(d) transmission des consignes des grandeurs de commande du procedd 
(Cqc) a des actionneurs, ou a des organes de reglage controlant des 
actionneurs, afin d'agir sur le deroulement du precede; 

dans laquelle 1'organe de prediction (OP) est base sur un module 
25 mathematique du precede, appete module direct (M), et est congu de telle 
man&re qu'on predit la masse M^r d'au moins un constituant (X) dans le 
reacteur (R) par F equation : 

M XR = LAG ^XRin • r X ■ r X > 

dans laquelle : 

" ^XRin est ^ e debit massique du constituant X entrant dans le reacteur R ; 
30 - t*x est le temps de sejour de x dans le rdacteur (constante de temps), qui 
vaut 

t x = / (E Fxdis) 

dans lequel : 

- M^R designe la derniere valeur calculee de la masse du 
constituant X presente dans le reacteur R ; 
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- £ Fxdis designe la somme de tous les debits rhassiques Fxdis avec 
lesquels le constituant X disparait du reacteur R, notamment par 
reaction et/ou par sortie du reacteur ; 
- la fonction y = LAG (u , r) est la solution de l'equation differentieUe 

u = x • ig * y 

5 calculee avec la valeur instantanee de u et de r, ainsi qu'avec la derniere 
valeur de y calcutee. 

L'avantage de cette methode est que l'equation differentielle 
susmentionnee est resolue par un calcul algebrique simple, par exemple par 
la formule suivante (T designant 1'intervalle de temps, generalement petit 
10 par rapport a t, separant les calculs successifs) ou par une formule 
equivalente a celle-ci : 

y(t-T) + u(t)-— 

y( t ) = ~ 



1 + 



T(t) 

Dans le cas ou les masses de plusieurs constituants sont evalu6es 
comme expose ci-dessus, la methode de I'invention est particulienement 
avantageuse dans la mesure ou ces masses peuvent etre calcul^es 

15 sequentiellement par de tels calculs algebriques simples, recalcules 
frequemment (T < < r en general). Par opposition, les methodes 
traditionnelles necessitent la resolution simultanee d'un systeme d'equations 
diff^rentielles, ce qui exige generalement une grande puissance de calcul et 
des algorithmes de calcul numerique (integration) sophistiques ; il en resulte 

20 que la dur£e de chaque iteration de calcul est elevee, et par consequent une 
regulation de ce type reagit mal a des variations rapides. 

Le precede de synthase regule peut servir a la synthase d'un 
compost monomerique ou polymerique ; de tres bons resultats ont 6t6 
obtenus pour la regulation de precedes de polymerisation. Le precede 

25 s'etend aussi au cas ou plusieurs produits utiles sont synthetises 

simultanement dans le meme precede. . Le precede peut etre continu ou 
discontinu (batch) ; la methode de regulation de r invention donne 
d'excellents resultats dans le cas de precedes continus. Le precede de 
synthese regule peut eventuellement ne consumer qu'une partie d'un 

30 precede plus vaste ? dont les autres parties sont regulees d'une maniere 

differente, ou ne sont pas regulees. Pour que la methode de regulation de 
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Tinvention soit applicable, il faut qu'au moins un reacteur puisse etre 
assimile a un reacteur parfaitement melange, c'est-a-dire a un reacteur dans 
Iequel les differentes grandeurs (temperature, concentrations des constituants 
presents, etc.) sont quasiment identiques en chaque point. D'autres 
5 reacteurs eventuels peuvent etre du type piston ("plug-flow") : on les 
modelise mathematiquement par des temps morts. La methode est 
egalement applicable a un procede se deroulant dans plusieurs reacteurs 
disposes en serie et/ou en parallele, pouvant produire des produits aux 
proprietes identiques ou differentes. 
10 Par "constituants", on entend designer Tensemble des substances 

presentes dans le reacteur et destinees a paiticiper h la synthase ou & la 
permettre : non seulement les reactifs de depart ainsi que le ou les produits 
synthetises, mais egalement la ou les £ventuelles substances ne subissant 
aucune transformation, telles que solvants, catalyseurs, etc. 
15 Outre un ou plusieurs reacteurs, Installation dans laquelle se 

derouie le procede peut eventuellement comprendre d'autres dispositifs 
classiques tels que detendeurs, strippers, condenseurs, secheurs, colonnes de 
distillation, etc. Generalement, ces dispositifs auxiliaires peuvent d'ailleurs 
etre egalement consideres comme des reacteurs (parfaitement melanges ou 
20 de type piston), meme si aucune reaction chimique ne s'y derouie. 

Dans le cas d'un procede de polymerisation, les "grandeurs liees aux 
proprietes du produit" peuvent par exemple etre choisies panni la masse 
moleculaire, l'indice de fluidite en fondu, la masse volumique standard, la 
teneur en comonomere lorsqu'un com onomfere est present, etc. 
25 Des exemples de "grandeurs liees au procede" sont notamment la 

temperature et la pression regnant dans le reacteur, failure de marche du 
procedd, les concentrations des differents reactifs dans le reacteur, etc. 
L'allure de marche designe la masse de produit synthetisee par unitd de 
temps, qui n'est cependant pas forcement egale au debit de produit 
30 synthetise sortant du reacteur : ainsi, a titre d'exemple, dans les phases de 
demarrage notamment, le debit massique de produit synthetise sortant du 
reacteur est generalement tr£s faible, voire nul, bien que la synthese ait 
commence, c'est-a-dire que ce debit sortant est alors inKrieur & Tallure de 
marche. En regime stable, par contre, on peut assimiler failure de marche 
35 a la masse de produit synthetise par unite de temps. 

Des exemples de "grandeurs de commande" sont les debits de 
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reactifs entrant dans le reacteur, la puissance fournie aux dispositifs de 
chauffage, etc. Ce sont ces grandeurs qui permettent d'agir sur le 
deroulement du precede ainsi que, generalement, sur les proprietes du 
produit synth&ise. 

5 La ou les consignes de la ou des grandeurs de commande sont 

transmises directement ou indirectement k des actionneurs classiques tels 
que notamment des vannes, des elements chauffants, etc. "Indirectement" 
signifie que les grandeurs de commande peuvent etre transmises par 
1' intermediate d'un ou plusieurs organes de regulation (controlant 

10 generalement une seule variable, par ex. des regulateurs PID) controlant le 
ou les actionneurs (regulation "locale"). 

Sur le plan materiel, l'organe de prediction et l'organe de controle 
sont generalement des dispositifs de calcul classiques permettant d'effectuer 
des calculs en fonction de leur cablage ou de leur programmation ; il peut 

15 notamment s'agir d'ordinateurs ou de systemes num£riques de controle- 

commande (SNCC). Un seul dispositif peut avantageusement combiner les 
fonctions de pi^diction et de contrdle. Le ou les dispositifs de calcul 
utilises sont de preference du type numerique, et foumissent periodiquement 
(par intermittence) les resultats de leurs calculs. Les intervalles de temps 

20 separant la fourniture de ces resultats peuvent varier dans le temps, et 

peuvent egalement differer selon le resultat concerne : il est clair que des 
grandeurs a variation rapide doivent etre recalculees plus frequemment que 
des grandeurs a variation lente. Des registres a decalage peuvent etre 
utilises pour simuler materiellement les temps morts. 

25 L'organe de prediction est base sur un modele mathematique direct 

du procede (M), dans lequel le reacteur (R) est assimile a un reacteur 
parfaitement melange ; un ou plusieurs delais purs (temps morts) peuvent 
eventuellement etre pris en compte pour representer d'eventuels reacteurs de 
type piston, d'eventuels delais de transport ou delais d'obtention de resultats 

30 de mesures, .... 

L'organe de controle est de preference base sur V inverse du modele 
direct utilise dans l'organe de prediction (modele inverse). 

En general, la somme E Fxdis de tous les debits massiques (Fxdis) 
avec lesquels le constituant X disparait du reacteur R comprend deux 

35 termes : 

* FR X qui designe le debit massique suivant lequel X est consomme par 
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une ou plusieurs reactions chimiques eventuelles ; 
* F Xout Q Ui d&i g ne Teventuel debit massique de X sonant du reacteur 
par soutirage en cours de reaction, dans le cas (habituel) ou X n'est pas 
entifcrement consomme par reaction dans ce reacteur ; ou encore, par 
5 exemple, par evaporation dans le cas d'un reacteur ouvert. 

L'interet de la methode est que les termes Fxdis sont en gendral 
proportionnels £ M XR > P ar exemple, on a en general 
F Xout = M XR 1 r R 

(t r designant le temps de sdjour du reacteur R) 
10 et Frx = R x . Mxr 

(R x designant la reactivite de X dans le reacteur R). 

Dans ce cas, P expression donnant r x se simplifie et devient : 

r x = 1 / (R x + l/r R ) 

Cette expression est independante de ce qui constitue une 

15 simplification extremement interessame. 

Un autre avantage de la methode reside dans le calcul periodique du 

temps de sejour r x . En effet, r x represente bien la dynamique du 

constituant consider^ dans le reacteur. Ceci permet notamment de suivre 

1'evolution de ce parametre, ce qui est important pour la comprehension de 
20 la dynamique du procede, et par consequent pour sa regulation. Au 

contraire, les methodes empiriques de type "boite noire" ne permettent pas 

d'acceder a ce parametre. 

Avantageusement, le calcul des predictions des grandeurs reglees 

(Pqr) peut additionnellement tenir compte d'une ou plusieurs mesures de 
25 grandeurs regimes (Mq R ), de grandeurs de commande (Mqq) et/ou d'autres 

grandeurs liees au ddroulement du procede (M^p). 

De meme, avantageusement, le calcul des consignes des grandeurs 

de commande du procede (Cqq) peut additionnellement tenir compte d'une 

ou plusieurs mesures de grandeurs reglees (Mq R ), de grandeurs de 
30 commande (Mq^) et/ou d'autres grandeurs liees au deroulement du procede 

(My^p), identiques ou diffdrentes de celles eventuellement prises en compte 

pour le calcul des predictions des grandeurs reglees (Pqr). 

Toutes les mesures dont il est question dans la presente description 

ne sont pas forcement des mesures directes, dans le sens ou Tune ou 
35 plusieurs d'entre elles peuvent eventuellement etre des mesures 

inferentielles. c'est-a-dire des valeurs obtenues par calcul a partir dune ou 
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plusieurs autres mesures directes. Ainsi, par exemple, Failure de marche 
de certains procedes de synthese exothermiques ne peut pas se mesurer 
directement, mais on peut en obtenir une mesure inferentielle par calcul, par 
exemple a partir de mesures (directes) du debit et des temperatures d'entree 
5 et de sortie du fluide de refroidissement. 

Dans le cas particulier de procedes de polymerisation, la ou les 
propriety du polymfcre intervenant dans la regulation sont de preference 
choisies parmi la masse volumique specifique (MVS) du polymere, les 
proprietes rheologiques du polymere a 1'etat fondu, et sa teneur en 

10 comonomere. En particulier, la ou les proprietes rheologiques intervenant 
dans la methode de regulation sont avantageusement l'indice de fluidite en 
fondu (melt index) du polymere et/ou une mesure de viscosity. 

Avantageusement, on evalue une ou plusieurs proprietes du polymfere 
en utilisant une technique choisie parmi la spectroscopic infrarouge proche 

15 (NIR), la spectroscopic infrarouge a transform ee de Fourier (FITR) et la 
resonance magnetique nucleaire (NMR). 

En particulier, on peut avantageusement evaluer une ou plusieurs 
proprietes du polymere en appliquant une relation de correlation preetablie 
aux resultats de mesures realises par spectroscopie infrarouge proche (NIR) 

20 & plusieurs longueurs d'onde predeterminees en fonction de ia nature du 
polymere, choisies entre 0,8 et 2,6 mm. 

Davantage de details concemant la realisation de telles mesures dans 
le cadre de la regulation de precedes de polymerisation peuvent etre trouves 
dans la demande de brevet EP 328826 (US 5155184). 

25 En vue de tenir compte d'ecarts eventuels entre les mesures et les 

predictions des grandeurs reglees, il peut etre utile de recourir a une 
correction. 

Un premier type de correction consiste en ce que la consigne d'au 
moins une grandeur reglee (C GR ) soit corrigee sur la base de l'ecait 

30 (avantageusement filtrS) entre la mesure (M GR ) et la prediction (Pqr) de 
cette grandeur reglee, de fa$on a ce que la regulation soit efficace (M GR = 
C GR ) meme en presence d'une erreur dans la prediction de cette grandeur 
nJglee. Cette technique est couramment designee par Texpression "internal 
model control" (IMC). 

35 Un second type de correction consiste en ce que le modele (M) du 

procede soit periodiquement adapte sur la base de Tecart (avantageusement 
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filtre) entre les predictions (Pqr) et les mesures (Mq R ) des grandeurs 
regl&s, de fagon a ce que, ici egalement, le modele du procede fournisse 
des predictions des grandeurs reglees (Pqr) aussi proches que possible 
(idealement, egales) des mesures de ces grandeurs (Mq R ), ce qui est 
5 indispensable a une regulation efficace. 

L'adaptation consiste h recalibrer le modele, c'est-a-dire a en 
recalculer un ou plusieurs parametres ; normalement, le nombre de 
parametres recalcules ne depasse pas le nombre des grandeurs reglees pour 
lesquelles on dispose a la fois d'une prediction et d'une mesure. Une 

10 resynchronisation (recalage dans le temps) de ces mesures est souvent 
souhaitable, surtout lorsqu'il s'agit de mesures de proprietes du produit 
synth&ise dont la duree d'obtention est longue. Ce second type de 
correction est plus avantageux dans la mesure bu il permet d 'adapter le 
modele egalement sur le plan de sa dynamique. 

IS L'adaptation conceme non seulement le modele direct du procede 

(organe de prediction), mais egalement le modele inverse (organe de 
controle). 

Selon une variante avantageuse, les mesures (Mq R ) des grandeurs 
reglees interviennent uniquement dans V adaptation eventuelle du modele du 
20 proc^d6, et n' interviennent pas directement dans le calcul des consignes des 
grandeurs de commande du procede (Cq^). 

C'est-k-dire que les mesures des grandeurs reglees n* interviennent 
pas dans la regulation proprement dite : Tavantage en est que la quality de 
la regulation n'est ainsi pas affectee par la lenteur eventuelle de revaluation 
25 des propriety du produit. 

Un autre aspect de l'invention concerne une methode de regulation 
telle que decrite ci-dessus, appliquee k un precede de polymerisation, 
comprenant Tune ou plusieurs des etapes supplementaires suivantes : 

calcul d'une consigne de temperature dans le reacteur en fonction d'une 
30 ou plusieurs consignes de proprietes du produit ; et transmission de 

cette consigne de temperature a un ou plusieurs actionneurs permettant 
de modifier la temperature dans le reacteur (eventuellement 
indirectement, e'est-a-dire par Tintermediaire d'un ou plusieurs organes 
de regulation, par ex. des regulateurs PID, controlant le ou les 
35 actionneurs) ; 

calcul d'un bilan thermique pour ie reacteur, sur base notamment de 



WO 96/41822 



VEP96/02408 



- 10- 



mesures de temperature; utilisation de ce bilan thermique de fagon a 
determiner la quantite de polymere synthetise par unite de temps (allure 
de marche) et/ou la productivity du catalyseur et/ou la concentration 
d'au moins un reactif dans le reacteur ; 
5 - calcul de la quantite de chaleur produite par la polymerisation, grace a 
un calcul de la quantite du ou des reactifs qui polymerisent ; 
determination par ce biais de la quantite la chaleur qu'il faut ajouter ou 
evacuer pour maintenir la temperature du reacteur ; utilisation du 
resultat dudit calcul (par exemple par feed-forward) pour ameliorer la 
10 regulation de temperature, de fa$on a respecter au mieux la consigne de 

temperature, notamment en cas de modifications de Failure de marche. 
Ces variantes sont basees sur le lien qui existe entre la quantite du 
ou des reactifs participant a la reaction et la quantite de chaleur produite ou 
absorbee par la reaction. 
15 Selon une variante avantageuse, la propriete Px R d'un constituant 

"x" dans le reacteur R, assimile a un reacteur parfaitement melange, est 
calculee comme suit : 

Px R = LAG (Pxjn , Mx R / Fxjn ) 
ou "Px w est une propriete d'un constituant "x", repondant de maniere 
20 substantielle a la loi de melange lineaire Pxi-j- 2 = w l - ^1 + w 2 ■ 

Px 2 , 

wj et W2 etant les proportions massiques de deux fractions 1 et 2 de 

propriete Pxj et Px 2 que Ton melange (wj + w 2 = 1) ; 
Px 1+2 est 13 P^op™te de x a sa sortie du reacteur, aprfes melange ; 
25 Pxin est la P'PP™ 1 ® du constituant M x" a son entree dans le reacteur R ; 
Mx R est la masse du constituant x dans le reacteur R ; 
Fxjjsj est le debit massique du constituant x entrant dans le reacteur R. 

Une transformation mathdmatique permet parfois de rendre lineaires 
(additives) certaines grandeurs qui ne le sont pas : par exemple, l'indice de 
30 fluidite en fondu d'un polymere (melt index) ne repond pas a une loi de 
melange lineaire, mais bien son logarithme ; le calcul de Px| +2 
susmentionne s'effecue des lors sur le logarithme de ce parametre. 

Selon une autre variante avantageuse, la methode de regulation de 
r invention comprend les etapes suivantes : 
35 - entree de consignes relatives a une ou plusieurs proprietes du produit a 
synthetiser, dans un algorithme maitre ; 
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- entree de la consigne d'allure de marche du procede, dans un algorithme 
esclave; 

- calcul des consignes de concentration des constituants dans le reacteur k 
Faide de Falgorithme maitre, en fonction notamment des consignes et des 

5 . mesures de proprietes du produit ainsi que de mesures ou de predictions 
des concentrations des diffiSrents constituants dans Ie rdacteur ; 

- transmission des consignes de concentration calculees par Falgorithme 
maitre comme grandeurs d'entree & Falgorithme esclave; 

- calcul de consignes de debit des constituants entrant dans le reacteur, a 
10 1'aide de Falgorithme esclave, en fonction notamment de la consigne 

d'allure de marche du procede, des consignes de concentration, et de 
mesures de ddbit des constituants entrant dans le reacteur, et 

- transmission des consignes de debit calculees a Faide de Falgorithme 
esclave a un ou plusieurs actionneurs (eventuellement indinectement, 

15 c'est-k-dire par Fintermediaire d'un ou plusieurs organes de regulation, 

par ex. des regulateurs PID, controlant le ou les actionneurs) afin de 
rtgler les debits des constituants entrant dans le reacteur, 
dans laquelle Falgorithme maitre et/ou Falgorithme esclave sont mis en 
oeuvre comme decrit ci-dessus, c'est-a-dire en utilisant la fonction LAG 
20 pour calculer la masse d'au moins un constituant dans le reacteur. 

Les algorithmes maitre et esclave sont dgalement mis en oeuvre au 
moyen d'un ou plusieurs dispositifs de calcul classiques. Selon une variante 
avantageuse, Fensemble des calculs (prediction, controle, etc.) de ces deux 
algorithmes sont effectues par un meme dispositif de calcul. 
25 Avantageusement, des mesures de temperature (par ex. temperature 

dans le reacteur, et/ou temperature d'entree et/ou de sortie d'un eventuel 
fluide de refroidissement) interviennent comme grandeurs d'entree 
supplementaires dans Forgane de prediction et/ou de controle. 

De preference, Falgorithme esclave prend £galement en compte des 
30 mesures de la composition de constituants presents dans ou sonant du 
reacteur. 

De maniere avantageuse, la methode de regulation comprend en 
outre une etape de calcul a Faide de Falgorithme esclave, en fonction des 
mesures de debits, de predictions de concentrations transmises a 
35 Falgorithme maitre pour calculer des predictions de proprietes servant 

comme grandeurs d'entree supplementaires dans le calcul des consignes de 



WO 96/41822 ^^ > CT/EP96/02408 



-12- 

concentrations. 

Les algorithmes maitre et esclave forment une regulation du type 
cascade. II est particulifcrement avantageux que l'algorithme maitre et/ou 
ralgorithme esclave soient adaptatifs, c'est-k-dire que certains de leurs 
5 parametres soient periodiquement recalcules (a intervalles reguliers ou 
imeguliers). Une telle adaptation permet notamment de garantir que le 
modele mathdmatique decrit le plus fidfelement possible le precede dans son 
etat courant, meme en cas de modification de certaines conditions 
operatoires (temperature, pression, allure de marche, etc.) et en cas de 
10 perturbations (empoisonnement du cataiyseur, ...). 

L'algorithme maitre realise le controle des proprietes du produit k 
I'aide d'un modele base sur des equations caracteristiques reliant les 
proprietes du produit aux concentrations des differents constituants dans le 
reacteur, ainsi qu'eventuellement a la temperature regnant dans le reacteur. 
15 L'algorithme esclave regule les concentrations d'un ou plusieurs constituants 
en agissant sur les debits d'alimentation d'un ou plusieurs constituants 
£ventuellement differents. 

L'interet de cette cascade "maitre-esclave" reside dans le fait que le 
modele maitre determine avec precision les concentrations de constituants 
20 necessaires a l'obtention des proprietes voulues pour le produit synthetise, et 
que le module esclave assure le respect des valeurs imposees par le maitre. 
Pilots par le maitre, Tesclave est par consequent capable : 
- d'amener rapidement les concentrations aux valeurs desirees par le maitre 
et de les y maintenir; 
25 -de controler efficacement Failure de marche du procede sans perturber 
les concentrations. 

Cette cascade maitre-esclave est particulierement efficace car aussi 
bien le maitre que V esclave prennent en compte la dynamique du procede 
grace a V utilisation de la fonction LAG dans les calculs. 
30 L'algorithme esclave peut en outre etre congu de maniere a fournir a 

l'algorithme maitre des predictions de concentrations fiables. De ces 
predictions ou de mesures de concentrations, ralgorithme maitre deduit des 
predictions de proprietes fiables du produit en train d'etre synthetise dans le 
rdacteur. En comparant ces predictions de proprietes aux consignes de 
35 proprietes, ralgorithme maitre peut, le cas echeant, intervenir et corriger 
les consignes de concentrations. Cette correction est possible avant meme 
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qu'un ecart entre une variable et sa consigne ne se presente. La prise en 
consideration de predictions de propriety obtenues a partir de predictions 
ou de mesures de concentrations permet de reduire considerablement les 
fluctuations temporelles des proprietes du produit synth£tis£ et il en resulte 
5 une meilieure Constance de la qualite du produit. 

Si les proprietes du produit a synthetiser sont d£pendantes de la 
temperature dans le ou les nSacteurs, il est preferable de prevoir une 
regulation de la temperature par 1'algorithme esciave. Ce dernier dresse ie 
bilan thermique de chaque reacteur et determine, au moyen du calcul de 
10 Failure de marche, la quantite de chaleur qu'il faut ajouter ou evacuer pour 
respecter des consignes de temperature calculees par 1'algorithme maitre. 
De ces resultats il ddduit des consignes d'entree pour des organes de 
regulation thermique de 1'appareillage de synthese. II sera apprecie que 
cette faqon de proceder permet d'intervenir sur les organes de regulation 
15 thermique de 1'appareillage de synthese avant meme que la temperature ne 
change. Des mesures de temperatures interviennent par ailleurs 
avantageusement comme grandeurs d'entree supptementaires dans 
1'algorithme esciave. 

L'algorithme maitre comprend avantageusement la structure 
20 suivante : 

- un organe de prediction, basd sur un modele direct du procede permettant 
de fournir une prediction des proprietes du produit synthetise en fonction de 
mesures et/ou de predictions des concentrations des constituants ; 

- un organe d' adaptation comparant les predictions de proprietes calculees 
25 par l'organe de prediction a des valeurs effectivement mesurees sur le 

produit synthetise et deduisant de cette comparaison des parametres 
d'adaptation, iesdits parametres d'adaptation intervenant comme grandeurs 
d'entree supplementaires dans ledit organe de prediction de Talgorithme 
maitre ; et 

30 - un organe de controle, base sur un module inverse du proced^, pour 
calculer, en fonction des consignes et des predictions de proprietes du 
produit a synthdtiser, des consignes de concentrations pour Falgorithme 
esciave, Iesdits parametres d'adaptation intervenant egalement comme 
grandeurs d'entree supplementaires dans ledit organe de controle. 

35 L'algorithme esciave comprend avantageusement la structure 

suivante : 
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- un organe de prediction, base sur un modele direct du procede permettant 
de fournir une prediction des concentrations d'un ou plusieurs des 
constituants sur base d'un bilan de matiere dans Ie reacteur ; 

- un organe d'adaptation, comparant les predictions de concentrations 

5 calculees par le modele direct a des mesures de concentration et deduisant 

de cette comparaison des parametres d' adaptation, lesdits paramfetres 
d * adaptation intervenant comme grandeurs d' entree supplementaires dans 
ledit organe de prediction de ralgorithme esclave ; et 

- un organe de controle, base sur un modele inverse du procede, pour 

10 calculer, en fonction de la consigne d'allure de marche, des consignes de 

concentration calculees par Forgane de controle de ralgorithme maitre et 
des predictions de concentration calculees par 1'organe de prediction de 
ralgorithme esclave, les consignes pour les debits entrant dans le 
reacteur, lesdits parametres d'adaptation intervenant comme grandeurs 

15 d'entree supplementaires dans ledit organe de controle de ralgorithme 

esclave. 

La dynamique du procede est avantageusement decrite et calculee au 
moyen de fonctions du type y = LAG (u , r) 

cette fonction etant la solution de l'equation differentielle 

dy 

u =T " dt +y 

20 ou les arguments u et t varient dans le temps. 1/ utilisation de cette 

fonction conformement aux thforemes 1 et 2 exposes ci-dessous permet de 
resoudre de maniere sequentielle les bilans de matiere utilises par 
ralgorithme esclave et de decrire la cinetique du procede par des equations 
caractdristiques simples dans ralgorithme maitre. La fonction LAG permet 

25 encore de reduire considerablement le volume des calculs necessaires et 
rend par consequent inutile l'utilisation d'ordinateurs rapides et puissants. 
De plus, cette fonction permet d'etablir de fagon particulierement simple les 
modeles directs et inverses du procede ou de certaines de ses parties. 
Les principals qualites de la regulation proposee peuvent etre 

30 resumees comme suit : 

- anticipation : la regulation commence a corriger les perturbations 
mesurees avant meme que leur effet ne se soit manifeste sur les mesures 
des proprietes (utilisation de predictions de concentrations, de predictions 
de proprietes et de predictions de temperatures dans les algorithmes); 
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- precision meme en presence de perturbations : le modele direct et le 
module inverse sont recalibres en permanence en utilisant les mesures des 
proprietes (adaptation); 

- validite etendue : I'algorithme garde sa validity durant les transitions 
d'allure de marche et de qualite, ainsi que lors des demarrages et des 
arrets (mise en equation de la dynamique du precede, utilisation de 
predictions pour les grandeurs dont les mesures impliquent des temps 
morts importants) ; 

- simplicite : Je developpement et la mise en oeuvre sont facility grace a 
une methode originale de mise en equation de la dynamique du procede 
(fonction LAG). 

Le procede de synthese a rgguler est ainsi modeiise sous une forme 
qui est generalement qualifiee de "modfele de connaissance" ("first principle 
model"), c'est-a-dire que son modele est 61abore a partir d'equations 
refletant le deroulement physico-chimique d6taill£ du procede. Une telle 
approche permet d'obtenir, au moyen d'un jeu d'equations relativement 
simple sur le plan mathematique, des rdsultats superieurs a ceux qui auraient 
6t6 obtenus au moyen d'un modele empirique de type boite noire, en 
foumissant notamment des parametres relies a des grandeurs reelles et une 
meilleure validite h l'exterieur de l'espace d' identification (extrapolation). 
La plupart des modeles empiriques utilisent des equations complexes, 
souvent d'ordre eleve si Ton desire obtenir une simulation correcte de la 
dynamique du precede, et dont les parametres (notamment les constantes de 
temps) doivent etre identifies pour un point de fonctionnement precis; le 
module n'est valable que dans le voisinage immgdiat de ce point de 
fonctionnement. Une telle approche peut difficilement etre generalisee k un 
grand nombre de points de fonctionnement dans le cas d'un procede de 
synthase chimique reel faisant intervenir de nombreuses grandeurs. 

Au contraire, selon la methode de regulation de V invention, on 
utilise un ensemble d'equations simples, purement statiques ; la dynamique 
du procede est simulee par des fonctions simples (cf. la fonction LAG ci- 
dessus). Avantageusement, les temps de sejour (constantes de temps des 
equations) peuvent etre recalcules aussi souvent qu'on le souhaite, ce qui ne 
pose aucun probleme etant donne la simplicite des equations. On obtient en 
fin de compte un ensemble d'equations extremement simples et faciles a 
resoudre en temps reel, meme a frequence elevee. 
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La methode de regulation proposee est avantageusement applicable a 
des precedes de synthese, notamment en continu, de poiymeres 
(polymerisation), et notamment a la polymerisation en continu d'olefines 
telles que par exemple 1'dthylene ou le propylene, aussi bien en phase 
5 liquide qu'en phase gazeuse. 

La presente invention concerne egalement un procede de synthase 
d'un ou plusieurs produits chimiques, regule au moyen d'une methode de 
regulation conforme a 1'invention. En particulier, de tres bons resultats ont 
ete obtenus pour la regulation d'un procede de synthese en continu de 

10 polyethylene par polymerisation d'ethylene dans au moins un reacteur, les 
reactifs comprenant l'ethylene, 1'hydrogene et/ou un comonomere optionnel, 
la reaction de polymerisation ayant lieu dans un solvant en presence d'un 
catalyseur et une partie du contenu du reacteur etant preleve en permanence 
ou par intermittence. Ce procede peut aussi bien se derouler en phase 

15 liquide qu'en phase gazeuse ; de preference, il se deroule en phase liquide 
(dans un solvant). 

La methode s' applique de maniere analogue a la synthese du 
polypropylene (le monomere de depart principal etant dans ce cas le 
propylene au lieu de l'dthylene), du propane pouvant egalement etre present 

20 si le procede se deroule en phase gazeuse. Pour le polypropylene, l'indice 
de fluidite en fondu est souvent designe par MFI au lieu de MI. 

L' invention concerne egalement un dispositif de regulation destine a 
mettre en oeuvre la methode de regulation de 1'invention, ainsi qu'une 
installation de synthese d'un ou plusieurs produits chimiques comprenant un 

25 tel dispositif de regulation. 

De maniere plus precise, Tinvention concerne egalement un 
dispositif de regulation d'un procede de synthese d'un produit chimique dans 
un appareillage de synthese comprenant au moins un reacteur, ledit 
dispositif comprenant : 

30 - au moins une unite de calcul; 

- des moyens pour entrer des consignes de proprietes du produit a synthe- 
tiser dans l'unite de calcul; 

- des moyens pour entrer des consignes d'allure de marche du produit a 
synthetiser dans l'unite de calcul; 

35 - des organes de mesure de debit des flux entrant dans le reacteur; 

- des organes de mesure de la composition des flux sortant du reacteur; 
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- des organes de commande de debits (actionneurs) pour regler les debits 
des flux entrant dans Ie reacteur; 

- des moyens de communication entre ladite unite de calcul, Iesdits organes 
de mesures de debit et Iesdits organes de regulation; 

S dans lequel : 

- la masse d'au moins un constituant est calculee par la fonction LAG 
comme exposd ci-dessus; 

- l'unit£ de calcul est apte k calculer a l'aide d'un aJgorithme maitre, en 
fonction des consignes de proprietes, des consignes de concentration de 

10 reactifs dans le reacteur; 

- r unite de calcul est apte a calculer a l'aide d'un algorithme esciave, en 
fonction des consignes de production et des consignes de concentration, 
de consignes de dSbit pour les flux entrant dans le reacteur, ces 
consignes de d£bit etant transmises comme consignes d'entree aux 

15 organes de regulation de debits; 

- les mesures realisees par les organes de mesure de debits interviennent 
comme grandeurs d'entree supplementaires dans ledit algorithme esciave, 
pour permettre a ce dernier de calculer, en fonction de ces mesures de 
debits, des predictions de concentrations, ces predictions de 

20 concentrations intervenant dans r algorithme maitre pour calculer des 

predictions de propriety utilisees comme grandeurs d'entree 
supplementaires dans le calcul des consignes de concentrations. 

L' invention conceme egalement un dispositif de regulation tel que 
decrit ci-dessus, dans lequel : 
25 - Tappareillage de synthese comprend en outre : 

des organes de regulation thermique apte a controler la temperature 

dans le reacteur, et 

des capteurs de temperatures; 

- Talgorithme maitre est apte a calculer, en fonction des consignes de 
30 propridtes, des consignes de temperatures pour le reacteur; 

- Falgorithme esciave est apte a : 

calculer un bilan thermique pour le reacteur, 
resoudre ce ou ces bilans thermiques de fagon a determiner la 
chaieur qu'U faut ajouter respectivement retrancher de la synthese 
35 pour respecter les consignes de temperature, et 

deduire de ce ou ces bilans thermiques des consignes d'entree pour 
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les organes de regulation thermique du reacteur, et 

recevoir comme grandeurs d'entree supplementaires les mesures 

effectuees par les capteurs de temperature. 

L'invention concerne egalement un dispositif tel que decrit ci-dessus, 
5 dans lequel des mesures effectuees par les capteurs de temperature 

interviennent comme grandeurs d'entree supplementaires dans ralgorithme 
maitre. 

L' invention concerne egalement un dispositif tel que decrit ci-dessus, 
comprenant : 

10 - au moins un analyseur apte k fournir des mesures des proprietes inter- 
venant dans I'algorithme maitre; et 

- des moyens pour entrer ces mesures de proprietes dans 1'unite de calcul; 
ladite unite de calcul comprenant : 

- un premier organe de prediction, base sur un premier modele direct du 
15 precede permettant la prediction des proprietes du produit synthetise, en 

fonction des predictions de concentrations calculees par Falgorithme 
esciave; 

- un premier organe d'adaptation, comparant les predictions de proprietes 
calculees par le premier organe de prediction aux valeurs effectivement 

20 mesurees sur le produit synthetise et deduisant de cette comparaison des 

parametres d'adaptation, intervenant comme grandeurs d'entree 
supplementaires dans ledit premier organe de prediction ; et 
. un premier organe de controle, base sur un premier modele inverse, pour 
calculer, en fonction des consignes et des predictions de proprietes, des 
25 consignes de concentrations pour l'algorithme esciave, lesdits parametres 

d'adaptation intervenant Egalement comme grandeurs d'entree 
supplementaires dans ledit premier organe de controle. 

L' invention concerne egalement un dispositif tel que decrit ci-dessus, 
comprenant en outre : 
30 - au moins un analyseur apte a fournir des mesures de concentration des 
reactifs; et 

- des moyens pour entrer ces mesures de concentration dans 1' unite de 
calcul; 

ladite unite de calcul comprenant : 
35 - un deuxifeme organe de prediction, base sur un deuxieme modele direct 
permettant la prediction des concentrations, en fonction du bilan de 
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matiere dans le reacteur; 

- un deuxieme organe d 'adaptation, comparant les predictions de 
concentrations calculees par le deuxieme organe de prediction aux 
mesures de concentration et deduisant de cette comparaison d'autres 

. parametres d'adaptation, intervenant comme grandeurs d'entrde 
supplementaires dans ledit deuxieme organe de prediction; et 

- un deuxieme organe de controle, bas6 sur un deuxieme module inverse, 
pour calculer, en fonction des consignes de production, des consignes de 
concentration calculees par Talgorithme maitre et des predictions de 
concentration du deuxieme organe de prediction, des consignes pour les 
debits entrant dans le reacteur, lesdits autres parametres d'adaptation 
intervenant comme grandeurs d'entree supplementaires dans ledit 
deuxieme organe de controle. 

Un mode de realisation concret de l'invention est illustre sur base 
d'un precede de synthase en continu de polyethylene (PE), en se referent 
aux Figures 1 a 10. Celles-ci montrent : 

Figure 1 un schema d'un circuit de fabrication de polyethylene; 

Figure 2 un schema simplifie de la structure d'une regulation 

avancee selon l'invention; 
Figure 3 un schema de principe de la regulation avancfe appliquee 

au circuit de fabrication de la Figure I; 
Figure 4 : un schema de principe d'un algorithme de controle 

adaptatif, tel qu' utilise dans le systeme de regulation 

avanc6e selon la Figure 2; 
Figure 5 un schema de la structure de 1'algorithme maitre dans le 

systfeme de regulation avancee selon la Figure 2; 
Figure 6 un schema de la structure de 1'algorithme esclave dans le 

systeme de regulation avancee selon la Figure 2; 
Figure 7 le schema general d'une m&hode de regulation conforme 

k l'invention; 

Fig. 8-10 le schema de variantes particulieres de la methode de 

l'invention. 

Sur la figure 7, on distingue tout d'abord le procede de synthese 
proprement dit (Pr), qui peut etre controle en imposant au moins une 
consigne d'une grandeur de commande (CGC) (par exemple un ou plusieurs 
debits de constituants entrant dans le reacteur. une temperature, etc.) k un 
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actionneur adequat (vanne, dispositif de chauffage ou de refroidissement, 
etc.). La regulation s'effectue au moyen d'un organe de controle (OC), 
base sur le module mathematique inverse du precede, et dont le role 
principal est de comparer la ou les consignes des grandeurs reglees (CGR) 
5 (par exemple une ou plusieurs proprietes du produit a synthetiser et/ou une 
ou plusieurs grandeurs liees au d^roulement du procede) a la ou aux 
predictions de ces grandeurs (PGR). La ou les predictions des grandeurs 
regimes (PGR) sont calculees par un organe de prediction (OP), base sur le 
module mathematique direct du procede, sur la base de mesures des 

10 grandeurs de commande (MGC). On notera qu'aucune mesure de 
proprietes) du produit synthetise n'intervient dans la regulation. 

La figure 8 reprgsente une variante de la methode de la figure 7, 
dans laquelle le module mathematique du procede est periodiquement adapte 
par un organe d'adaptation (OA), sur la base de 1'dcart (avantageusement 

15 filtre, ou traite numeriquement) entre les predictions (PGR) et les mesures 
(MGR) des grandeurs reglees. Une resynchronisation (recalage dans le 
temps) de ces mesures et de ces predictions est souvent necessaire, par 
exemple lorsqu'il s'agit de mesures de proprietes du produit synthetis£ dont 
la duree d'obtention est longue. L'organe d'adaptation (OA) transmet le 

20 resultat de ses calculs, c'est-a-dire ses instructions d'adaptation, a l'organe 
de prediction (pour adapter le modele direct du procede) et a l'organe de 
controle (pour adapter le modele inverse du procede). On notera que la ou 
les mesures de propriete(s) du produit synth&ise ne sont pris(es) en compte 
que dans le processus d'adaptation, qui se deroule generalement a une 

25 frequence beaucoup plus faible que le processus normal de controle. La 
lenteur eventuelle de ces mesures n'a done aucun effet direct sur la qualite 
de la regulation. 

La figure 9 represente une autre variante de Finvention, dans 
laquelle une ou plusieurs mesures de grandeurs reglees (M GR ) sont prises 

30 en compte par l'organe de controle (OC), et une ou plusieurs mesures de 
grandeurs reglees (M QR ') (eventuellement differentes) sont prises en 
compte par l'organe de prediction (OP). De meme, une ou plusieurs 
mesures de grandeurs de commande (Mqq) peuvent en outre etre prises en 
compte par l'organe de controle (OC). 

35 II va de soi qu'on pourrait creer une autre variante de la methode de 

1' invention en combinant les variantes des figures 8 et 9, c*est-a-dire a la 
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fois en utilisant un organe <T adaptation et en prenant en compte une ou 
plusieurs mesures de grandeurs de commande dans F organe de controle, 
et/ou une ou plusieurs grandeurs de commande dans Torgane de prediction 
et/ou dans F organe de controle. 

Dans la figure 10, le modele mathematique du procdde n'est pas 
adapte au sens propre du terme, mais on utilise Pecan (avantageusement 
filtre) entre les mesures et Ies predictions des grandeurs reglees pour 
corriger les consignes des grandeurs r£gl6es (C GR ). En F occurrence, cette 
correction a 6t6 ici representee comme une simple difference : on soustrait a 
chacune des consignes de grandeurs reglees un terme correctif calcule par 
F organe d'adaptation OA (qui n'est en fait ici qu'un organe de correction), 
ce qui fournit des consignes corrigees C GR \ transmises a F organe de 
controle OC. II va de soi que dans certains cas, la correction peut 
impliquer des operations plus complexes qu'une soustraction, par exemple 
une division (dans ce cas-ci, on peut cependant se ramener a une 
soustraction en considerant les logarithm es des grandeurs considerees). 
Cette methode est communement appelde Internal Model Control (IMC). 

Si Ton se refere a la Figure 1, qui represente sch£matiquement un 
circuit de synthese en continu de polyethylene (PE), la polymerisation de 
Fethylene se fait dans un reacteur boucle 10, en suspension dans un solvant 
adequat, tel que par exemple 1'hexane. Le processus est continu, c'est-a- 
dire que les reactifs sont injectes en continu, et qu'une partie du contenu du 
reacteur 10 est continuellement soutiree. Une pompe de circulation (non 
representee) assure l'homogeneit£ du contenu du reacteur 10. 

Les reactifs introduits dans le reacteur sont Tethylene "Et", 
Thydrogene "Hy" et le butene "Bt" (cf. reference 11). Un catalyseur est 
egalement injecte en continu. II importe de bien maitriser les concentrations 
des reactifs dans le reacteur, car les propriet£s de la resine PE sont 
principalement determinees par les rapports de concentrations Hy/Et et 
Bt/Et. 

La temperature de polymerisation dans le reacteur est un parametre 
additionnel qui influence les proprietes de la resine PE. Vu que la reaction 
de polymerisation est fortement exothermique, la temperature du reacteur 
doit etre regulfe en ayant recours a un circuit de refroidissement 12. 

Le reacteur 10 en fonctionnement contient done du solvant, du poly- 
mere et des reactifs n'ayant pas encore reagi et du catalyseur. Son contenu 
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est soutire en permanence a travers la conduite de soutirage 14. Ce contenu 
soutire entre dans un stripper STP 16, qui sdpare Ie polymere PE et les 
fluides (solvant et reactifs). Ces fluides sont vaporises par injection de 
vapeur d'eau et evacues dans un condenseur CD 18. Dans ce dernier le 
5 solvant est de nouveau condense avant d'etre recycle. Les reactifs qui sont 
plus lagers sont separes du solvant et sont aussi recycles. Un chromato- 
graphe en phase gazeuse GC (20) place a la sortie du condenseur 18 permet 
de determiner les concentrations Hy/Et, Bt/Et des reactifs. 

Le polymere evacu6 du stripper 16 est concentre dans une 

10 centrifugeuse CFG 22, puis sech£ dans un secheur a lit fluide SHLF 24, 

avant d'etre envoye vers le finishing pour y etre granule. Des echantillons 
sont prelevds a la sortie du secheur 24 pour y mesurer les proprietes de la 
resine : cristallinite (mesuree par la masse voliimique specifique "MVS") et 
proprietes rheologiques (indice de fluidite en fondu ("melt index (MI)" ou 

15 "melt flow index (MFI)") et viscosite en fondu V 2 M , mesuree dans un 
ecoulement capillaire sous un cisaillement de i00 s~*). 

La dynamique de ce processus de synthase de PE est lente et 
complexe : 

- Le reacteur boucle 10 se comporte comme un reacteur parfaitement 
20 melange. Des lors, toute modification du debit d'alimentation d'un des 

reactifs ne se traduira que progressivement sur la concentration de ce 
reactif dans le reacteur. En effet, le nouveau debit doit se melanger au 
volume tout entier du reacteur 10 pour 1'amener a la nouvelle concen- 
tration d'6quilibre. 

25 - La mesure des concentrations des reactifs se fait par un chromatographe 
en phase gazeuse 20; il s'agit d'un appareil discontinu, qui procfede par 
etapes successives : prelevement d'un echantillon de gaz, analyse, puis 
elaboration des resultats. II y a done un temps mort (de 5 a 15 min.) 
entre les changements de concentration et leur mesure. 

30 - Les proprietes du polymere fabrique a chaque instant dependent principa- 
lement des concentrations des reactifs. Toute modification de ces 
concentrations affecte done instantanement les proprietes du polymere 
fabrique. Par contre, les proprietes moyennes dans le reacteur ne se 
modifient que progressivement, car Ie polymere fraichement produit doit 

35 se melanger au polymere deja present dans le reacteur 10 (temps de 

sejour ; ± 2 h). 
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- Lorsque le poly mere est soutirg du reacteur 10, il subit a nouveau une 
serie de melanges dans les differents appareils (STP, CFG el SHLF) 
destines a le secher (temps de sejour : ± 2 h). Des echantillons de 
polymere sont alors preleves et analyses par le laboratoire de 1'usine. 
5 . Les resultats de ces analyses ne seront done communiques & la fabrication 
qu'apres un nouveau temps mort, qui peut etre important (± 2 h). 
Modelisation du procedg h Taide d'une fonction LAG 

Selon la methode de l'invention. la modelisation dynamique d'un 
precede de synthese en continu est realisee en ayant recours & des hypo- 
10 theses de melanges parfaits et de delais purs. Les melanges parfaits sont 
mis en equation au moyen d'une fonction bien connue des ingenieurs, la 
fonction "LAG", ou filtre passe-bas (du ler ordre) : cette fonction est 
lineaire et aisement programmable. 
Elle est definie comme suit : y = LAG (u , r) 
15 (dit "LAG de u durant r") comme etant la solution de V equation 
differentieLle 

u = x • + y 

dont les arguments u et r varient dans le temps. 

Cette Equation peut se resoudre numeriquement (meme en temps 
reel) par une equation algebrique du ler ordre, ayant pour arguments les 
20 variables suivantes : 

- la periode d'echantillonnage "T" (ou temps ecoule depuis la derniere 
iteration) 

- le temps de sejour (ou "constante" de temps) V a Finstant "t" 

- la grandeur d'etat "y" & V instant precedent "t-T" 

25 - la grandeur de commande "u" a Tinstant present "t" 

(u et r representent en fait les valeurs mesurees ou calculees a T instant "t" 
sur les grandeurs "u" et "r", que Ton suppose avoir ete constantes dans tout 
1'intervaUe precedent "T"). 

De preference, T est petit par rapport a t (par exemple au moins 10 
30 fois plus petit), de fa^on a augmenter la precision du calcul. 

La solution de Tequation susmentionnee peut se calculer par exemple 
par la formule suivante : 

y(t) = y(t-T) . e " T/r (t) + u(t ) . ( l . e " T/r ( t > ) 
ou encore, plus simplement (approximativement) : 
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La modelisation du precede a l'aide de la fonction LAG se base sur 
Ies thdoitmes suivants : 

Soit un reacteur parfaitement melange (CSTR), de volume V R . 
Divers constituants (reactifs ou inertes) l'alimentent, parnii lesquels le 
5 reactif "x" (debit entrant Fx^) ayant la propriete "Px^" a Tentree. On 
mesure egalement le debit de sortie "F OUT " (soutirage). 
Theoreme 1 : Application de la methode LAG au calcul d'un bilan 
massique : 

A chaque instant, la masse M Mx R " d'un constituant "x" dans un 
10 reacteur parfaitement melange (CSTR) est egale au LAG du produit du debit 
massique entrant "Fxjj^" par un temps "r x ", durant ce meme temps "r x " : 
Mx R = LAG (Fx m . r x , r x ) (en kg) 

Le temps M r x " est le "temps de sejour de x"; il vaut la masse du constituant 
Mx R , divisee par la somme des debits massiques "sonants" (quantite cons- 
15 ommee par reaction "FRx", debit soitant du reacteur w Fx OUT ", etc.). 

r x = Mx R / (Fx QUT + FRx + .*.) (en h ) 

Ce theoreme donne ainsi une methode exacte pour le calcul dynamique 
(meme en temps rdel) des concentrations dans un reacteur parfaitement 
melange. En effet, la concentration "Cx R " du constituant "x" exprimee en 
20 kg/m 3 vaut, si "V R " est le volume du reacteur : 

Cx R = Mx R / V R (en kg/m 3 ) 

Connaissant le debit volumique total "FV OUT " sortant du reacteur, 
on definit le temps de sejour du reacteur "r R " : 

T R - V R 1 ^OUT 
25 Dfcs lors, le d<Sbit massique M Fx OUT " du constituant "x" en sortie du 

reacteur (soutirage) vaut : 

Fx OUT = Mx R 1 r R 

On remarquera que si "x" est un inerte (ne subissant pas de reaction 
et ne sortant du reacteur que par soutirage), on a : 
30 r x = r R 

Par ailleurs, dans le cas frequent ou la vitesse de reaction de "x" est 
proportionnelle a sa concentration Cx R , avec un facteur de proportionnalite 
R x , on a : FRx = R x . Mx R et des lors 
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r x = 1 / (R X + l/r R ) 
Thdoreme 2 : Application de la methode LAG au calcul d'une propriete 
d'un melange : 

Soit une propriete "Px" d'un constituant, repondant k la loi de 
5 melange lineaire : 

Pxj + 2 = w l • P*i + w 2 * ***2 
ou wj et W2 sont les fractions massiques de propriete Pxj et Px 2 . (avec Wj 
4- w 2 = 1). 

A chaque instant la propriete Px R dans un reacteur parfaitement 
10 melange (CSTR) est £gale au LAG de la propriete a r entree Pxjj^, durant 
un temps de sejour qui vaut le rapport de la masse Mx R du constituant dans 
le reacteur divise par le debit massique entrant (et/ou apparaissant par 
reaction) Fxjj^ : 

Px R = LAG ( Px m , Mx R / Fx^ ) 
15 II est ainsi possible de tenir compte de la dynamique du proced6, et 

de recalculer en permanence ses constantes de temps. 

Comme expose precedemment, les "propriet6s" dont il est question 

ici peuvent dans certains cas avoir subi une transformation mathematique les 

rendant lineaires (par exemple, le logarithme de Findice de fluidity en fondu 
20 d'un polymere peut etre considere comme repondant h une loi de melange 

lineaire). 

Principe de la regulation : 

Le modele du procede de synthase etant etabli, il faut un algorithme 
qui calcule les parametres necessaires a la regulation de ce proc6d£. 

25 La Figure 2 montre un schema general simplifie du type de 

regulation avanc6e (Advanced Process Control ou "APC") adopte pour le 
procede de polymerisation d^crit plus haut. On voit que ce systeme de 
regulation comprend une cascade de deux algorithmes, cene cascade pilotant 
notamment les regulateurs PID des debits d' alimentation en reactifs. 

30 Les deux algorithmes en cascade, appeles algorithme maitre et algo- 

rithme esclave, sont tous deux des algorithmes dynamiques adaptatifs bases 
sur des modeles issus de la connaissance du processus (par opposition aux 
modeles empiriques), bases notamment sur les bilans de matiere et la 
cinetique du procede regule. lis utilisent la fonction LAG predefinie. 

35 Sur la Figure 3, qui montre le principe du systeme de regulation 

dans le contexte du procede de polymerisation decrit plus haut, on voit 



WO 96/41822 



7EP96/02408 



-26- 



que : 
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15 



20 



25 



30 



- I'algorithme maitre est base sur les equations caracteristiques des cata- 
lyseurs, c'est-a-dire des equations fournissant les proprietes du PE en 
fonction de la temperature de polymerisation et des concentrations des 
reactifs dans le reacteur; il founiit a I'algorithme esclave les consignes 
des concentrations en reactifs pour satisfaire les consignes pour les 
proprietes du PE; 

- 1'algorithme esclave est base sur un bilan matifere et la cinetique chimique 
des reactions; il fournit aux regulateurs des debits d' alimentation les 
consignes de debits de reactifs necessaires pour satisfaire les consignes de 
concentrations imposees par I'algorithme maitre et la consigne d'allure de 
marche du procede. De preference, il calcule egalement un terme 
anticipatif ("feed-forward") pour la consigne de temperature, ameliorant 
la regulation de temperature lors des changements d'allure de marche.^ 

Ce type de regulation n'est d'une precision parfaite que si le modele 
est parfait et prend en compte toute les perturbations possibles. Cela n'est 
en general pas le cas. C'est pourquoi en general (voir Figure 4) on adapte 
en continu le module direct (et le module inverse) en comparant les 
predictions aux mesures des proprietes. Cette "adaptation" du modele 
permet de maintenir sa precision en presence de perturbations non 
modelisees, et ainsi d'obtenir une regulation plus precise en toutes 
circonstances. 
L'alporithme esclave 

La Figure 6 illustre le principe de 1'algorithme esclave : 

1. un organe de prediction utilisant un module direct du procede predit, h 
partir des mesures des debits alimentant le r6acteur, les concentrations 
en reactifs et en polymere; 

2. un organe d' adaptation compare les concentrations en ethylene (Et), 
hydrogene (Hy) et butene (Bt) mesuree par un analyseur (chroma- 
tographe a gaz) aux valeurs predites par le modele direct, de fagon a 
determiner trois parametres d' adaptation : 

- l'activite specifique du catalyseur pour l'ethylene "AEt", en kg/h de 
polyethylene par kg de catalyseur et par kg/m 3 d'ethylene 

- reireur de gain sur la mesure de debit d'hydrdgene "KfHy" 

- la purete de Talimentation en butene "KfBt" 

3. T organe de controle calcule, a panir des consignes de concentration 
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calculees par Falgorithme maitre et de la consigne d* allure de marche, 
les consignes pour Ies debits d'alimentation des reactifs ; ces consignes 
de debit sont composees d'un feed-forward base sur le modele inverse, 
et d'un feed-back pnoportionnel a l'ecart entre le modfele direct et les 
5 consignes de concentration. 

Pour comprendre les calculs effectues par ralgorithme esclave, dans 
le cas d'un precede de synthese de polyethylene, il faut d'abord savoir qu'il 
est generalement admis que la vitesse de polymerisation "VitP^ 0 est 
proportionnelle : 

10 - a la concentration en Et non polymerise cEt R (en kg/m 3 ), 

- a la concentration en catalyseur actif dans le reacteur cCA R (en kg/m 3 ), 
et 

- a un facteur de proportionnalite, fonction (difficilement quantifiable) de 
la temperature, des concentrations en Hy, en Bt et en co-catalyseur, de la 

15 presence de poisons, etc. Ce terme est appele "activite catalytique" pour 

F&hylene Agj. Sauf dysfonctionnement majeur (poison, etc.), il varie 
raisonnablement peu en cours de campagne. 
vitP Et = A Et • cCA R . cEt R (kg/m 3 .h) 

La quantite d'&hylene "FpEt" qui polymerise dans le rdacteur a 
20 chaque unite de temps ("ddbit de polymerisation") vaut done, si V R est le 
volume du reacteur : 

FpEt « VitP^ . V R (kg/h) 

= Agj . cCA R . MEt R 
ou MEt R est la masse d' ethylene en solution dans le reacteur (en kg) 
25 On sait par ailleurs que la vitesse d' incorporation de FHy est environ 

100 fois plus lente que FEt, et 10 fois plus lente pour le Bt. On en deduit : 
FpHy = Agj . cCA R . MHy R / 100 
FpBt » A m . cCA R . MBt R / 10 

ou MHy R est la masse d'hydrogene en solution dans le reacteur (en kg) et 



30 MBt R est la masse de butene en solution dans le reacteur (en kg) 
Le modele esclave utilise maintenant les mesures suivantes : 

FEtj N = debit d' alimentation en ethylene (monomere) (kg/h) 

FSv IN = d ^ bit d ' alimentation en solvant (hexane) (kg/h) 

FCAjjq- debit d 'alimentation en catalyseur (kg/h) 

35 FHy IN = debit d'alimentation en hydrogene (agent de transfert) (kg/h) 

FBt IN = debit d' alimentation en butane (co-monomere) (kg/h) 
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II utilise egalement les parametres cT adaptation suivants : 
A El = "activite catalytique" pour F ethylene 

KfHy = erreur de gain sur la mesure de Falimentation en hydrogfcne 
KfBt = purete de 1' alimentation en butfcne. 
5 Les calculs suivants sont effectues de fagon sequentielle, dans Tordre 

suivant, avec une frequence elevee (l'intervalle de temps separant chaque 
iteration etant faible par rapport au plus petit des temps de sejour r x ). 

Le volume du reacteur etant constant, Ie debit volumique sortant 
equivaut au debit volumique entrant (fluides incompressibles). Le debit 
10 volumique sortant "FVout" P eut d ^ s lore etre calcule comme la somme 

des debits massiques entrants divises par la masse volumique qu'ils ont dans 
le reacteur : 

FV OUT = FSv^/650 + FEtjN/950 + FBt^/600 (m 3 /h) 
(ou les masses volumiques sont les suivantes : 650 kg/m 3 pour le solvant, 

15 950 kg/m 3 pour le polyethylene, 600 kg/m 3 pour le butene). On fait ici 
l'hypothfese que tout Fethylfene se transforme instantanement en poly- 
ethylene, et on neglige le debit d'hydrogfene et de catalyseur (quelques kg). 

Le solvant est chimiquement inerte et ne sort du reacteur que par 
soutinige du reacteur. On calcule sa masse "MSv R " dans le reacteur en 

20 utilisant le theoreme 1 : 

Tr _ v R / FV OUT (h) (temps de sejour dans le reacteur) 

MSv R = LAG (FSv^ . r R , t r ) (kg) 

Le catalyseur se desactive avec une constante de temps "kd" ; la 
masse "MCA R " de catalyseur actif dans le reacteur est calculee comme 

25 suit : 

r CA = 1 / (l/r R + kd) (h) 
MCA R = LAG (FCA m . t ca , t ca ) (kg) 
et done la concentration "cCA R " de catalyseur actif dans ie reacteur vaut : 
cCA R = MCA R / V R (kg/m 3 ) 
30 L'&hylene "sort" du reacteur dans le debit de soutirage ainsi que 

dans la reaction de polymerisation. Sa masse dans le reacteur "MEt R " est 
calculee comme suit : 

r Et = 1 / (l/r R + A £t . cCA R ) 00 
MEt R = LAG (FEtjjsj . r Et , r Et ) (kg) 
35 D'une fagon similaire, on calcule la masse "brute" (non calibree) 

d'hydrogene dans le reacteur "MHy RAW " : 
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r Hy = 1 / (l/r R + Agj . cCA R / 100) (h) 

MHyRAW = LAG (WlN • r Hy - 7 Hy ) ^) 
La masse "MHy R ", corrigee pour prendre en compte l'erreur de 

gain sur la mesure de F alimentation en hydrogene, vaut : 

5 MHy R = KfHy . MH yRAW (kg) 

Le rapport Hy:Et dans le reacteur vaut done : 

HyEt R = MHy R / MEt R 

D'une fagon similaire, on calcule la masse de butene "brute" 

10 r Bt = 1 / (l/r R + A m . cCA R / 10) (h) 
MBtj^w = LAG (FBtflvj . T Rt , r Rt ) (kg) 

La masse M MBt R ", corrigee pour prendre en compte la purete de 
F alimentation en butfene, vaut : 

MBt R = KfBt . MBtj^w (kg) 
15 Le rapport Bt:Et dans le reacteur vaut done : 

Bt£t R = MBt R / MEt R 

On a vu que le debit de polymerisation "FpEt" (allure de marche 
instantanee) vaut : 

FpEt = A £t . cCA R . MEt R (kg/h) 
20 Le polymfcre 6tant inerte et ne subissant pas de reaction, sa masse 

dans le reacteur M MPE R " vaut : 

MPE R = LAG ( FpEt . r R , r R ) (kg) 

Le debit de polymere sortant du reacteur "FPEqux" vaut des lors : 
FPEqut - MPE R / r R (kg/h) 

25 Adaptation du modele esclave : 

Le bloc d'adaptation utilise un analyseur (par exemple un chroma- 
tographe a gaz "GC") pour obtenir les mesures des concentrations dans le 
reacteur en Ethylene "cEt GC ", en hydrogene w cHy GC " et en butene 
"cBt GC " (exprimees par exemple en kg/m 3 ). Ces mesures sont comparees 

30 aux valeurs predites par le module direct, afin de determiner les trois 
parametres d'adaptation suivants : 

- F activity speScifique du catalyseur pour F ethylene "AEt", en kg/h de 
polyethylene par kg de catalyseur et par kg/m 3 d'ethylene 

- Ferreur de gain sur la mesure de debit d'hydrogene "KfHy" 
35 - la purete de Falimentation en butene "KfBt" 

Le chromatogTaphe a gaz fournit des mesures echantillonnees avec 
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un retard d'environ 6 minutes. 

Le calcul de 1'activite specifique pour 1'ethylene "AEt" decoule des 
egalites suivantes : 

cEt GC = cEt R (kg/m 3 ) 
5 = MEt R / V R 

= LAG (FEt IN . r Et , r R ) / V R 
= LAG (FE^ / (l/r R + A £t . cCA R ) , r m ) / V R 
=* Agj" = ( FEtjjsj / LEAD ( cEt QC . V R , r Et ) - 1/t r ) / cCA R 

La mesure cEt GG etant 6chantillonnee et bruitee, on ne souhaite pas 
10 la trouver dans un LEAD; des lors, on utilise plutot : 
A et ' = LAG (A^" ■ T Et> 

» LAG (FEt IN /cCA R , T m ) I (cEt GC .V R ) - 1 / LAG (r R .cCA R , r a ) 
On tient compte du retard de 6 minutes sur la mesure de cELqq en 
introduisant deux LAGs en serie de 3 minutes chacun sur les valeurs du 
15 module, et on obtient la formule finale : 

A et - LAG ( FEtjj^ / cCA R , r Et , 3/60 , 3/60 ) / (cEt GC . V R ) 
- 1 / LAG ( r R . cCA R , r a , 3/60 , 3/60 ) 

Le calcul du gain du debit d'hydrogene "KfHy" decoule des egalites 
suivantes : 

20 cHy GC - cHy R (kg/m 3 ) 
= MHy R / V R 
- KfHy . MH yRAW / V R 
KfHy= cHy GC . V R / MH yRAW 

On tient compte du retard de 6 minutes sur la mesure de cHy G £ en 
25 introduisant deux LAGs en serie de 3 minutes chacun sur la valeur du 
modele, et on obtient la formule finale : 

=> KfHy = cHy QC . V R / LAG ( MHy RAW , 3/60, 3/60 ) 

D'une fagon similaire a KfHy, on calcule le parametre de correction 
pour la purete du butene "KfBt" : 
30 KfHy = cBt GC . V R / LAG ( MBtj^w , 3/60, 3/60 ) 
L'algorithme de regulation a pour entrees : 
- les consignes de concentration calculees par l'algorithme maitre; plus 
precisement les consignes pour les rapports de concentration cHy R /cEt R 
"HyEt SP " et cBt R /cEt R "BtEt sp " (en kg/kg) 
35 - la consigne pour Failure de marche du procede FpEtgp, fixee par 
Toperateur 
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- la consigne pour la concentration en Ethylene cEt sp , fixee par 
Toperateur; 

- les concentrations calculees par le modele. 

II calcule les consignes pour les debits d'alimemation des reactifs 
5 ^^SP' FCA SP» ^ySP et ^^SP* Dif ^ ren ts algorithmes peuvent etre 
utilises, parmi lesquels le MBPC (Model Based Predictive Control). lis 
peuvent en general se decomposer en un feed-forward base sur le module 
inverse, et un feed-back proportionnel k Pecan entre le modele direct et les 
consignes de concentration. 
10 Regulation de P alimentation en Ethylene : 

Feed-forward : valeur pour maintenir la concentration actuelle, basee sur 
T inversion de la valeur stationnaire de V equation suivante : 

MEt R = LAG (FEt^ . , ) * (kg) 
FEtpp = MEt R / r a 
15 Feed-back : proportionnel a Tecart entre la consigne cEt sp , et le modfele 
FEtpg = 5 . (cEt sp . V R - MEt R ) 
Consigne : 

FEt sp = FEtpp + FEtpg 
Regulation de Falimentation en catalvseur : 
20 Feed-forward : valeur pour maintenir la concentration actuelle, basee 

sur 1'inversion de la valeur stationnaire de r equation suivante : 

MCA R - LAG (FCA IN . r QA , r CA ) (kg) 
FCApp = MCA R / r CA 
Feed-back : proportionnel a Tdcart entre la consigne FpEt sp , et le module, 
25 selon la formule suivante : 

FpEt = A m . MEt R . MCA R / V R 
=* FCApg = 5 . (FpEt sp / (Apj . MEt R / V R ) - MCA R ) 
Consigne : 

FCA sp = FCApp + FCApg 
30 Regulation de Talimentation en hvdrogene : 

Feed-forward : valeur pour maintenir la concentration actuelle, basee 
sur Tinversion de la valeur stationnaire de Tequation suivante : 

MHy R = LAG (FH yiN . r Hy , r Hy ) (kg) 
FHypp = MH yRAW / r Hy 
35 Feed-back : proportionnel h I'ecart entre la consigne de rapport HyEt Sp , et 
le modele. 
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FHypB = 5 . (HyEt Sp . MEt R - MHy R ) 
Consigne : 

FHy sp = FHypp + FH ypB 
Regulation de V alimentation en butene (similaire a l'hydrogene) : 

FBtpp = MB tRAW / r Bt 

FBtpg = 5 . (BtEtgp . MEt R - MBt R ) 

FBt Sp = FBtpp + FBtpg 

Les equations qui precedent resument les Equations de ralgorithme 
esclave. Elles sont executees toutes les 10 secondes par le systfcme nume- 
rique de controle et de commande (SNCC). 
L'algorithme maitre 

La Figure 5 illustre le principe de ralgorithme maitre : 

1 . son organe de prediction (base sur un modele direct) predit les 
proprietes principales du polymere (MI et MVS); il utilise pour cela la 
mesure de la temperature de polymerisation, les predictions pour les 
concentrations dans le reacteur, fournies par le modele esclave, et les 
temps de sejour du PE dans les difterents appaieils ; 

2. son organe d'adaptation compare les mesures de MI et de MVS 
effectutes (k la sortie du secheur egalement), soit par le laboratoire de 
mesure toute les 2 h, soit par un analyseur en continu, aux valeurs 
predites par le modele direct, de fagon a determiner les 2 parametres 
d'adaptation, qui sont des parametres correctifs, multiplicatif pour le MI 
et additif pour la MVS ; 

3. son organe de controle (base sur un module inverse) calcule, a partir des 
consignes de MI et de MVS foumie par Toperateur, les consignes pour 
les concentrations dans le reacteur (rapports Hy:Et et Bt:Et). Comme 
pour ralgorithme esclave, ce calcul est compose d'un feed-forward base 
sur le module direct, et d'un feed-back proportionnei a 1'ecart entre le 
modele direct et les consignes de l'operateur. 

Pour un catalyseur donne, les proprietes des resines en regime stable 
sont des fonctions de la temperature de polymerisation et des concentrations 
des reactifs. Parmi les differentes Equations statiques renseignees dans la 
litterature, on a choisi les equations suivantes : 
log(MI) = aQ + a! . T° + a 2 . log (Hy/Et) + a 3 . Bt/Et 
MVS = b Q + bj . T° + \>2 . (Bt/Et) b3 + b 4 . log(MI) 

Les parametres aQ a a 3 et bg a b 4 sont obtenus par identification en 
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regime stable, pour plusieurs resines fabriquees avec le meme catalyseur. 

Par ailleurs, les differents appareils rencontres par le polyethylene 
jusqu'au moment ou ses proprietes sont mesurees (reacteur, stripper, 
centrifugeuse puis secheur), peuvent tous etre assimiles en premiere appro- 
5 ximation a des reacteurs parfaitement melanges. 

L'algorithme maitre utilise comme entrees les mesures suivantes : 
T R = temperature au reacteur (°C) 
Vstp = volume liquide dans le stripper (obtenu par mesure de niveau) (m J ) 
^MES = mesure de MI (melt-index) 
10 MVSj^jes = mesure de MVS (Masse Volumique Sp&ifique) 
ainsi que les calculs suivants realises par ralgorithme esclave : 
FpEt = production instantanee de polymere (allure de marche) (kg/h) 
FPE OUT = debit de PE sortant du reacteur (kg/h) 
JMQPE R = masse de PE dans le reacteur (kg) 
15 * ilyEt R = rapport Hy sur Et dans le reacteur (kg/kg) 
BtEt R = rapport Bt sur Et dans le reacteur (kg/kg) 

Les valeurs instantanees brutes (avant adaptation) de la MVS et du 
logarithme du MI ("IMF) sont calcul^es par : 
MVS^ = b 0 + bj . T R + b 2 . (BtEt R ) b3 + b 4 . lMI^g 
20 IMIjns = aQ + a t . T R -h a 2 . log (HyEt R ) + a 3 . BtEt R 

En utilisant le theoreme 2, on calcule les proprietes moyennes brutes & 
la sortie du reacteur : 
IMIr = LAG (IMI^s, MPE R / FpEt) 
MVSr = LAG (MVS INS , MPE R / FpEt) 
25 En effet : 

- les proprietes 1MI et MVS repondent assez bien a une loi de melanges 
lineaire 

- le reacteur boucle peut etre assimile a un reacteur parfaitement melange 

- le debit massique de PE "entrant" (apparaissant) dans le reacteur est 
30 effectivement FpEt, la quantity de PE qui polymerise a chaque instant 

(allure de marche). 

Proprietes brutes a la mesure : Sachant qu'on a dans le stripper 
environ 500 kg de PE par m 3 , et si Ton prend Fhypothfcse que le stripper 
est un reacteur parfaitement melange, on calcule les proprietes brutes en 
35 sortie du stripper comme suit : 

LMIstp = LAG (IMIr, 500 . Vstp / FPE OUT ) 
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MVSstp = LAG (MVSr, 500 . Vstp / FPE QUT ) 

Le temps de sejour dans la centrifugeuse etant tres court, on peut 
l'ignorer. 

Le s&heur est a lit fluide; il contient en permanence environ 
5 1. 400 kg de PE. On peut faire Thypothfese que le niveau dans le stripper 

change peu, et que le d6bit qui en sort egale celui qui y rentre. Des lors, le 

debit de PE entrant dans le secheur vaut FPE OUT . On a alors en sortie du 

secheur, a l'endroit ou l'echantillon est preleve pour la mesure de propriete, 

les valeurs brutes suivantes : 
10 IMIsh = LAG (lMIstp, 1400 / FPE OUT ) 

MVSsh = LAG (MVSstp, 1400 / FPE OUT ) 

Les proprietes apres adaptation sont obtenues en faisant intervenir les 

param&tres d'adaptation kMI (parametre multiplicatif) et kMVS (parametre 

additif); les proprietes apres adaptation sortie reacteur, stripper et secheur 
15 valent done : 

MIr c = kMI . lO 1 ^ 

MVSr c = kMVS + MVSr 

MIstp c = kMI . 10^^ 

MVSstp c = kMVS + MVSstp 
20 MIsh c = kMI . 10>MIsh 

MVSsh c = kMVS + MVSsh 

Adaptation de Talgorithme maitre : 

Les mesures de proprietes prennent un certain temps a etre realisees 

(± 5 min. si analyseur en ligne, ± 1 h si effectuees par le laboratoire). 
25 Pour pouvoir calculer les param^tres d'adaptation, il faut done 

resynchroniser (recaler dans le temps) les predictions brutes du modfele avec 

les mesures. Cela peut par exemple se faire au moyen d'un registre a 

decalage (appele ici "fonction DELAY") : 

1MI DEL = DELAY (IMIsh, tj^j) 
30 MVS DEL = DELAY (MVSsh, t mvs ) 

avec : tj^j et r MVS = ± 5 min. ou ± 1 h, suivant que la mesure est 

realisee par un analyseur continu ou par le laboratoire. 

A chaque nouvelle mesure de MI ou de MVS, on recalcule le 

parametre d'adaptation brut kMI' ou kMVS', en comparant la valeur 
35 modele brute resynchronisee a la valeur mesuree : 

kMI' = log (MImes) - IMlDEL 
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kMVS' - MVS MES - MVS D£L 

Ces valeurs brutes sont filtrees afin d'attenuer les reactions 
precipitees que les perturbations (bruit) eventuelles des mesures infligeraient 
au precede : 
5 kMI = LAG ( kMF, ± 1 h) 

kMVS = LAG ( kMVS\ ± 1 h) 
Organe de controle : 

L'organe de controle a pour consignes les valeur MI S p et MVSgp, 
entrees par Toperateur. II calcule les consignes pour les rapports des 
10 concentrations dans le reacteur HyEt sp et BtEt sp necessaire pour obtenir 
rapidement les proprietes voulues MIgp et MVS<jp. Ce calcul se fait en 2 
etapes : 

1. Torgane de controle calcule, a partir des consignes Mlgp et MVSgp 
fournies par 1'operateur, et des valeurs aprfcs adaptation du MI et de la 

15 MVS dans les difterents appareils, les consignes Mligp et MVSigp pour 

la production instantanee. Ces consignes instantanees sont composees 
d'un feed-forward, et d'un feed-back proportionnel k F6cart entre le 
modele direct et les consignes de l'op^rateur. 

2. les consignes pour les rapports de concentration HyEtgp et BtEt SP sont 
20 alors calculees en inversant V equation statique utilisee plus haut pour le 

calcul de la valeur instantanee du MI et de la MVS. 

Consignes pour les proprietes instantanees : On compare les 
proprietes sortie secheur aux consignes de proprietes, pour determiner les 
consignes voulues pour les proprietes a la sortie du stripper (on neglige la 
25 centrifugeuse) : 

MIstp sp = 10 < »og(MISP) + 0.1 . (log(MISP) - log(MIshC))) 

MVSstp SP = MVS sp + 0.1 . (MVS Sp - MVSsh c ) 

De meme, a partir de l'ecart entre ces consignes sortie stripper et 

des valeurs calibrees au stripper, on calcule les consignes voulues en sortie 

30 du reacteur : 

Mir = 10 < log(MIstpSP) + 0.5 . (log(MIstpSP) - log(MIstpC))) 

MVSr Sp = MVSstp sp + 0.5 . (MVSstp Sp - MVSstp c ) 

Enfln, a partir de Tecart entre ces consignes sortie reacteur et des 

valeurs calibrees correspondantes, on calcule les consignes voulues pour la 

35 production instantanee : 

mi = 10 < l°g(MIrSP) + 2 . (log(MIrSP) - log(MIrC))) 
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MVSi SP = MVSr sp + 2 . (MVSr sp - MVSr c ) 

Consignes pour les rapports de concentrations : Les consignes pour 
les rapports de concentration HyEtgp et BtEtgp s'obtiennent en inversant 
l'equation statique utilisee plus haut pour le calcul de la valeur instantanee 
5 du MI et de la MVS, en substituant aux termes MI et MVS les consignes 
voulues pour la production instantanee et en appliquant le parametre 
d* adaptation. 

Partant de : 

log (Mli sp / kMI) = 2lq + a! . T R + a 2 . log (HyEt sp ) + a 3 . BtEt R 
10 MVSi sp - kMVS = b Q + bj . T R 4- b 2 . (BtEt sp ) b3 + b 4 . lMI^g 
on obtient : 

a 2 - log (HyEt sp ) = log (MIi sp / kMI) - (a^ + aj . T R + a 3 . BtEt R ) 
ce qui donne : 

HyEt sp = io ( ( lo S (MKSP / kMI) - aO - al . TR - a3 . BtEtR ) / a2 ) 
15 et : 

b 2 . (BtEt Sp ) b3 = MVSi sp - kMVS - ( b Q + bj . T R + b 4 . IMI^g ) 
ce qui donne : 

BtEt sp = ( ( MVSi Sp - kMVS - b Q - bj . T R - b 4 . lML ms ) / b 2 ) 1/b3 
Les equations qui precedent resument les equations de V algorithm e 
20 maitre. Elles sont executees toutes les 30 secondes par le SNCC. 

Avec ce procede, il est possible de piloter la polymerisation avec une 

grande precision. En particulier : 

- les proprietes controlees (MI et MVS) sont maintenues au plus pres des 
valeurs voulues, avec une dispersion minime 

25 - les changements de qualite (et done des proprietes MI et MVS) sont 
realises avec rapidity et precision 

- les demarrages et arrets de la polymerisation, ainsi que les changements 
d 'allure de marche du procede, sont realises d'une fagon acc£l£ree, tout 
en maintenant le MI et la MVS fort proches des valeurs voulues. 

30 Bien que la methode de regulation selon 1' invention ait ete presentee 

a 1'aide d'un procede de synthese de polyethylene par polymerisation en 
continu d'ethylene, il est entendu que cette methode de regulation sera de 
fagon generate efficace pour d'autres procedes de synthese, et en particulier 
pour les procedes presentant Tune ou plusieurs des caracteristiques 

35 suivantes : 

- une regulation multivariate est necessaire parce que plusieurs variables 
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influent sur I'ensemble des proprietes a regler: 

- la dynamique du processus est lente : melanges en serie. temps morts 
importants; 

- les mesures des propridtes sont echantillonnees avec une faible frequence, 
et/ou bruitees; 

- Ia regulation doit etre dynamique, c'est-a-dire valable quelle que soit 
Failure de marche du proced£, ainsi que durant les transitions d*allure de 
marche et de qualite (proprietes) du produit a synthetiser; 

- il est interessant d'estimer certaines variables non directement mesurees. 
Pour pouvoir etre facilement mise en oeuvre avec les techniques 

presentees, il suffit que : 

- les equations statiques du processus soient connues (souvent, elles le sont 
au moins dans une certaine mesure, sinon le procede ne pounait etre 

m maitris< * 

15 - la dynamique du processus puisse etre approchee par des melanges 
parfaits et des temps morts; 

- les mesures necessaires soient disponibles et d'une quality suffisante (en 
particulier, les debits des reactifs et les debits traversant les reservoirs 
concernes). 

20 L'utilisation particulifere de la fonction LAG d^crite ci-dessus, 

notamment dans les theoremes 1 et 2, peut naturellement etre Vendue a des 
methodes de regulation basees sur une structure autre que celle qui a ete 
exposee, qui comprend un algorithme maitre et un algorithme esclave 
distincts. EUe peut par exemple etre appliquee dans une methode de 

25 regulation ne comportant qu'un algorithme unique. 
Exemple? : 

8 essais de synthase de polyethylene (PE) de 4 types differents 
(ddfinis par leur MI, MVS, ...) ont ete realises, respectivement en utilisant 
une mdthode de regulation classique et en utilisant la methode de 
30 Tinvention. Le tableau ci-dessous resume les constatations qui ont 6t6 faites 
sur la base de nombreuses mesures de 1'indice de fluidite en fondu (melt 
index) des 8 polymeres obtehus. Cpk designe Tindice de capability centre 
du procede. 
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Type 
dePE 


Regulation classique 


Regulation selon l'invention 


Ecart-type 


Cpk 


Ecart-type 


Cpk 


1 


0.127 


0.909 


0.059 


2.202 


2 




0.6i 




2.0 


3 




0.48 




1.88 


4 




0.64 




1.09 



On constate que 1'indice de capability Cpk est plus que double grace 
a r utilisation de la methode de l'invention, ce qui indique que les proprietes 
sont environ deux fois moins dispersees et/ou mieux centrees par rapport 
aux consignes imposees. 
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Tableau des abreviations utilisees 
aj parametres pour F equation statique du MI (I = 0 a 3) 

bj parametres pour F equation statique de la MVS (I = 0 & 4) 

Agj activite catalytique pour F ethylene (m^.kg'^h" 1 ) 

5 cx GC concentration de M x" obtenue par la mesurc de Fanalyseur (kg/m 3 ) 

cx R concentration de "x" dans le reacteur (kg/m 3 ) 

cx§p consigne pour la concentration de "x" dans le reacteur (kg/m 3 ) 

Fpx debit massique de polymerisation de "x" (allure de marche) (kg/h) 

FV OUT debit volumique sonant du reacteur (m 3 /h) 

10 Fxjj^j debit massique de "x" entrant (kg/h) 

Fx OUT d * Jt massi q ue de " x " sonant (kg/h) 

kd constante de desactivation du catalyseur (1/h) 

KfBt parametre correctif (d'adaptation) pour le butene 

KfHy parametre correctif (d'adaptation) pour Fhydrogene 

15 kMI parametre correctif (d'adaptation) pour le MI 

kMVS parametre correctif (d'adaptation) pour la MVS 

LAG ( , ) fonction filtre passe-bas du ler ordne 

M^MES mesure du MI (Melt Index) (indice de fluidite en fondu) 

MIy MI (Melt Index) brut (non calibnS) dans w y" 

20 MIy c MI calibre (avec adaptation) dans M y" 

MtysP consigne de MI (calibree) pour "y H 

MVS MES mesure de la MVS (Masse Volumique Standard) 

MVSy MVS (Masse Volumique Standard) brute (non calibrS) dans "y M 

MVSy c MVS calibree (avec adaptation) dans M y" 

25 MVSy sp consigne de MVS (calibree) pour M y" 

Mx RAW masse de "x" brute (non calibre) dans le reacteur (kg) 

Mx Y masse de "x" calibree (avec adaptation) dans "y" (kg) 
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reactivite de X dans le reacteur 
V Y Volume de "y" (m 3 ) 

r R temps de sejour dans le reacteur (h) 

temps de sejour pour "x M dans le reacteur (h) 
5 "x" peut representer les constituants suivants : 
Bt Butene 
CA Catalyseur 
Et Ethylene 
Hy Hydrogene 
10 Sv Solvant 

"y" peut representer les appareils suivants : 
r Reacteur de polymerisation 

stp Stripper 
sh Secheur (a lit fluide) 
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Leeende des figures 

10 reacteur de polymerisation 

11 alimentation en reactifs (matieres premiferes), catalyseur, solvant 

12 circuit de refroidissement 
5 14 conduite de soutirage 

16 stripper 
18 condenseur 

20 chromatographe en phase gazeuse 
22 centrifugeuse 
10 24 secheur a lit fluide 

25 solvant et reactifs a recycler 

26 polyethylene 

27 reactifs a recycler 

28 solvant a recycler 

15 30 consignes de propriete du polymere 

31 consigne d' allure de marc he du precede 

32 algorithme maitre 

33 algorithme esclave 

34 consignes de concentration 
20 35 consignes de debits entrants 

36 regulation (PID) des debits 

37 mesures 

38 regulation de temperature 

39 feed-forward de temperature 

25 40 dynamique de la polymerisation : cin&ique chimique et bilan matiere 

41 simulations de l'ailure de marche et des rapports Hy/Et et Bt/Et 

42 mesures de temperature, de debits et de concentrations 

43 procede regule 

44 analyse d'un echantillon du produit synthetise par le precede 
30 51 mesure de grandeurs liees au deroulement du procede 

52 mesure de proprietes du polymere 

53 modele direct ; prediction des proprietes a la mesure 

54 comparaison : correction du module (adaptation) 

55 algorithme de regulation, base siir le modele inverse (feed-forward + 
35 feed-back) 
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56 consignes pour les grandeurs liees au precede 

57 mesures et consignes de debits entrants 

58 mesure de la temperature et prediction des concentrations 

59 modele direct : equations des proprietes en fonction des concentrtations 
5 60 prediction des proprietes 

61 consignes pour ies concentrations dans le reacteur 

62 mesure des debits entrants 

63 mesure des concentrations dans le reacteur 

64 modele direct ; prediction des concentrations sur base du bilan matiere 
10 65 predictions de concentration 

66 comparaison ; calcul des parametres d' adaptation 
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REVENDIC ATIONS 

1 - Mgthode de regulation d'un proc&te de synthese d'au moins un 
produit chimique dans un appareillage comprenant au moins un reacteur (R) 
pouvant etre assimile k un reacteur parfaitement melange, dans laquelle une 
5 • ou plusieurs grandeurs de commande (GC) permettent d'agir sur le 

deroulement du precede en vue qu'une ou plusieurs grandeurs liees aux 
proprietes du produit et/ou au deroulement du precede, appelees grandeurs 
regldes (GR), soient dgales k des consignes corespondantes (Cq R ), ladite 
methode comprenant les etapes suivantes : 
10 (a) entree de consignes concernant les grandeurs reglees (Cq R ) ; 

(b) calcul, au moyen d'un organe de prediction (OP), de predictions des 
grandeurs reglees (Pqr)» sur * a base de mesures des grandeurs de 
commande du precede (Hqq) ; 

(c) utilisation d'un organe de controle (OC) pour calculer des consignes 
15 des grandeurs de commande du precede (Cqq), sur base des 

consignes (C GR ) et des predictions (Pqr) des grandeurs reglees ; 

(d) transmission des consignes des grandeurs de commande du proc&te 
(Cqq) k des actionneurs, ou a des organes de reglage controlant des 
actionneurs, afin d'agir sur le deroulement du precede ; 

20 dans laquelle I'organe de prediction (OP) est base sur un module 

mathematique du precede, appele modele direct (M), et est congu de telle 
maniere qu'on pr&iit la masse d'au moins un constituant (X) dans le 

reacteur (R) par V equation : 

M XR = lAG ^XRin • r X ■ r X > 

25 dans laquelle : 

- FxRin est * e d€b\t massique du constituant X entrant dans le reacteur R ; 

- est le temps de sejour de x dans le reacteur, qui vaut 

T X = M XR 1 < E Fxdis ) 

dans lequel : 

30 - M XR d^signe * a demiere valeur calculee de la masse du constituant 

X presente dans le reacteur R ; 

E Fxdis designe la somme de tous les debits massiques Fxdis avec 
lesquels le constituant X disparait du reacteur R, notamment par 
reaction et/ou par sortie du reacteur ; 
35 - la fonction y = LAG (u , r) est la solution de T equation differentieUe 
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dy 

calculee avec la valeur instantanee de u et de r, ainsi qu'avec la derniere 
valeur de y calculee. 

2 - Methode de regulation selon la revendication 1 , dans laquelle la 
consigne d'au moins une grandeur regime (Cqr) est corrigee sur la base de 

5 l'ecart entrc la mesure (Mq R ) et la prediction (Pgr) de cette grandeur 

reglee, de fagon a ce que la regulation soit efficace meme en presence d'une 
erreur dans la prediction de cette grandeur regime (Pgr). 

3 - M6thode de regulation selon la revendication 1 , dans laquelle le 
module (M) du precede est periodiquement adapte sur la base de 1'ecart 

10 entre Ies predictions (Pgr) et les mesures (Mqr) des grandeurs reglees, de 
fagon a ce que le modele du procede fournisse des predictions des grandeurs 
rdglees (Pgr) aussi proches que possible des mesures de ces grandeurs 
(M GR ). 

4 - Methode de regulation selon la revendication 3, dans laquelle les 
15 mesures (M G r) des grandeurs reglees interviennent uniquement dans 

r adaptation eventuelle du modele du proc&le, et n'interviennent pas 
directement dans le calcul des consignes des grandeurs de commande du 
procede (Cqq). 

5 - Methode de regulation selon Tune des revendications 1 a 4, 
20 appliquee k un procede de polymerisation, comprenant Tune ou plusieurs 

des etapes supplementaires suivantes : 

- calcul d'une consigne de temperature dans le reacteur en fonction d'une 
ou plusieurs consignes de proprietes du produit ; et transmission de cette 
consigne de temperature a un ou plusieurs actionneurs permettant de 

25 modifier la temperature dans le reacteur ; 

- calcul d'un bilan thermique pour le reacteur, sur base notamment de 
mesures de temperature; utilisation de ce bilan thermique de fagon a 
determiner la quantite de polymere synthetise par unite de temps et/ou la 
productivity du catalyseur et/ou la concentration d*au moins un reactif 

30 dans le reacteur ; 

- calcul de la quantite de chaleur produite par la polymerisation, grace a un 
calcul de la quantite du ou des reactifs qui polymerisent ; determination 
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par ce biais de la quantite la chaJeur qu'il faut ajouter ou evacuer pour 
maintenir la temperature du reacteur ; utilisation du resultat dudit calcul 
pour ameliorer la regulation de temperature, de fagon a respecter au 
mieux la consigne de temperature, notamment en cas de modifications de 
5 . Failure de marche. 

6 - Methode de regulation selon Tune des revendications 
precedentes, dans Iaquelle la propria Px R d'un constituant "x" dans le 
reacteur R, assimile a un reacteur parfaitement melange, est calculde comme 
suit : 

10 Px R = LAG (Pxjn , Mx R / Fx^ ) 

ou "Px" est une propriete d'un constituant "x", repondant de maniere 
substantielle a la loi de melange lineaire Pxj+j = w l ■ P*i + w 2 • Px 2 • 
W| et W2 etant les proportions massiques de deux fractions 1 et 2 de 
propria Px| et Px 2 que Ton melange ; 

15 Pxj + 2 ^ l a propri&d de x a sa sortie du reacteur, apr&s melange ; 

Pxjj^j est la proprtete du constituant "x" a son entree dans le rdacteur R ; 

Mx R est la masse du constituant x dans le reacteur R ; 

Fx IN ^ ' e debit massique du constituant x entrant dans le reacteur R. 

7 - Mdthode de regulation selon Tune des revendications precedentes, 
20 comprenant les etapes suivantes : 

- entree de consignes relatives a une ou plusieurs proprietes du produit a 
synthetiser, dans un algorithme maitre ; 

- entrde de la consigne d 'allure de marche du precede, dans un algorithme 
esclave; 

25 - calcul de consignes de concentration des constituants dans le reacteur k 

l'aide de Talgorithme maitre, en fonction notamment des consignes et des 
mesures de proprietes du produit ainsi que de mesures. ou de predictions 
des concentrations des diffSrents constituants dans le reacteur ; 

- transmission des consignes de concentration calculees par 1'algorithme 
30 maitre comme grandeurs d'entree a 1'algorithme esclave; 

- calcul des consignes de debit des constituants entrant dans le reacteur, a 
l'aide de l'algorithme esclave, en fonction notamment de la consigne 
d'allure de marche du precede, des consignes de concentration, et de 
mesures de d^bit des constituants entrant dans le reacteur, et 

35 - transmission des consignes de debit calculees a l'aide de r algorithme 
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esclave a un ou plusieurs actionneurs afin de regler les debits des 

constituants entrant dans le reacteur, 
dans laquelle l'algorithme maitre et/ou l'algorithme esclave sont mis en 
oeuvre conformement a Tune des revendications prec&ientes. 

5 8 - Methode de regulation selon la revendication 7, caract£risee en 

ce que Falgorithme maitre comprend : 

- un organe de prediction, base sur un module direct du procede permettant 
de fournir une prediction des proprietes du produit synthdtise en fonction 
de mesures et/ou de predictions des concentrations des constituants ; 

10 - un organe d'adaptation comparant les predictions de proprietes calcul6es 
par T organe de prediction a des valeurs effectivement mesurees sur le 
produit synthetise et deduisant de cette comparaison des parametres 
d'adaptation, lesdits parametres d'adaptation intervenant comme 
grandeurs d'entree supplementaires dans ledit organe de prediction de 

15 l'algorithme maitre ; et 

- un organe de controle, base sur un module inverse du procede, pour 
calculer, en fonction des consignes et des predictions de proprietes du 
produit a synthetiser, des consignes de concentrations pour l'algorithme 
esclave, lesdits parametres d'adaptation intervenant egalement comme 

20 grandeurs d'entree supplementaires dans ledit organe de controle. 

9 - Methode de regulation selon Tune quelconque des revendications 
7 ou 8, dans laquelle l'algorithme esclave comprend : 

- un organe de prediction, base sur un modele direct du procede permettant 
de fournir une prediction des concentrations d'un ou plusieurs des 

25 constituants sur base d'un bilan de matiere dans le reacteur ; 

- un organe d'adaptation, comparant les predictions de concentrations 
calculdes par le modele direct k des mesures de concentration et 
deduisant de cette comparaison des parametres d'adaptation, lesdits 
parametres d'adaptation intervenant comme grandeurs d' entree 

30 supplementaires dans ledit organe de prediction de Talgorithme esclave ; 

et 

- un organe de controle, base sur un modfele inverse du procede, pour 
calculer, en fonction de la consigne d'allure de marche, des consignes de 
concentration calculees par !' organe de controle de Talgorithme maitre et 

35 des predictions de concentration calculees par 1' organe de prediction de 
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Talgorithme esclave, les consignes pour les debits entrant dans Ie 
reacteur, lesdits parametres d'adaptation intervenant com me grandeurs 
d'entree supplementaires dans Jedit organe de controle de Talgorithme 
esclave. 

5 10 - Methode de regulation selon Tune des revendications 

precedentes, appliquee a la regulation de la synthese en continu de 
polyethylene par polymerisation d 'Ethylene dans au moins un reacteur, les 
reactifs comprenant ['ethylene, 1'hydrogene et/ou un comonomere optionnel, 
la reaction de polymerisation ayant lieu en presence d'un catalyseur et une 
10 partie du contenu du reacteur etant preleve en permanence ou par 
intermittence. 

1 1 - M&hode de regulation selon les revendications 9 et 10, dans 
laquelle Torgane d'adaptation de Talgorithme esclave compare les mesures 
des concentrations en ethylene (Et), hydrogene (Hy) et/ou comonomfere 

15 optionnel (Bt) aux valeurs predites par Toj^gane de prediction de 

Talgorithme esclave, de fagon k determiner au moins Tun des parametres 
d'adaptation suivants : 

a) 1'activite specifique du catalyseur pour l'&hylfene "AEt", en kg/h de 
polyethylene par kg de catalyseur et par kg/m 3 d'ethylene ; 
20 b) 1'erreur de gain sur la mesure de debit d'hydrogene "KfHy" ; 
c) la purete de 1'alimentation en comonomere "KfBt". 

12 - Methode de regulation selon Tune des revendications 1 a 9, 
appliquee a la regulation de la synthase en continu de polypropylene par 
polymerisation de propylene dans au moins un reacteur, les reactifs 

25 comprenant le propylene, 1'hydrogene et/ou un comonomere optionnel, la 
reaction de polymerisation ayant lieu en presence d'un catalyseur et une 
partie du contenu du reacteur etant preleve en permanence ou par 
intermittence. 

13 - Methode de regulation selon les revendications 9 et 12, dans 
30 laquelle 1'organe d'adaptation de Talgorithme esclave compare les mesures 

des concentrations en propylene (Pe), hydrogene (Hy) et/ou comonomere 
optionnel (Et) aux valeurs predites par 1'organe de prediction de 
Talgorithme esclave, de fagon a determiner au moins Tun des parametres 
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cT adaptation suivants : 

a) Tactivite specifique du catalyseur pour le propylene "APe", en kg/h de 
polypropylene par kg de catalyseur et par kg/m 3 de propylene ; 

b) Temeur de gain sur la mesure de debit d'hydrogene "KfHy" ; 
5 c) la puretd de I'alimentation en comonomere "KfEt". 

14 - Methode de regulation selon la revendication 8, appliquee k un 
precede de polymerisation, dans laquelle : 

- on mesure periodiquement l'indice de fluidite en fondu (MI) et/ou la 
masse volumique standard (MVS) du polymere et/ou sa teneur en 

10 comonomere; 

- Torgane de prediction de l'algorithme maitre calcule des predictions 
brutes de MI et MVS en fonction de la temperature au reacteur, des 
concentrations dans le reacteur et des temps de sejour dans les differents 
appareils du circuit de polymerisation ; 

15 - periodiquement, Forgane d'adaptation de Talgorithme maitre : 

resynchronise les predictions brutes de MI et MVS, en tenant compte 
du temps ecoule entre la prise des mesures de MI et MVS et 
l'obtention du resultat des mesures, et compare les predictions brutes 
de MI et MVS resynchronisdes aux mesures de MI et MVS, 

20 - calcule un parametre d'adaptation multiplicatif kMI applique a la 

prediction brute du MI pour obtenir une prediction calibree du MI, 
et 

calcule un parametre d'adaptation additif kMVS applique a la 
prediction brute de MVS pour obtenir une prediction calibree de 
25 MVS. 

15 - Methode de regulation selon Tune des revendications 
precedentes, appliquee a un procede de polymerisation, dans laquelle on 
evalue une ou plusieurs proprietes du polymfere en utilisant une technique 
choisie parmi la spectroscopic infrarouge proche (NIR), la spectroscopie 
30 infrarouge a transformee de Fourier (FITR) et la resonance magnetique 
nucleaire (NMR). 
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16 - Methode de regulation selon Tune des revendications 
precedentes, appliquee a un procede de polymerisation, dans laquelle on 
evalue une ou plusieurs proprietes du polymere en appliquant une relation 
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de correlation pre&ablie aux rgsultats de mesures realisees par spectroscopic 
infrarouge proche (NIR) a plusieurs longueurs d'onde predetermin&s en 
fonction de la nature du polymere, choisies entre 0,8 et 2,6 mm. 

17 - Procede de synthase d'un produit chimique dans un appareillage 
5 • comprenant au moins un reacteur pouvant etre assimile a un reacteur 

parfaitement melange, reguld au moyen de la mdthode de regulation selon 
Tune des revendications precedentes. 

18 - Dispositif de regulation d'un procede de synthese d'un produit 
chimique dans un appareillage de synthese comprenant au moins un reacteur 

10 pouvant etre assimile a un reacteur parfaitement melange, dans lequel la 
m&hode de regulation est conforme a Tune des revendications 1 & 16. 
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